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Abstract 

Weak limits as the density tends to infinity of classical ground states of integrable pair potentials 
are shown to minimize the mean-field energy functional. By studying the latter we derive global 
properties of high-density ground state configurations in bounded domains and in infinite space. Our 
main result is a theorem stating that for interactions having a strictly positive Fourier transform the 
distribution of particles tends to be uniform as the density increases, while high-density ground states 
show some pattern if the Fourier transform is partially negative. The latter confirms the conclusion 
of earlier studies by Vlasov (1945), Kirzhnits and Nepomnyashchii (1971), and Likos et al. (2007). 
Other results include the proof that there is no Bravais lattice among high-density ground states of 
interactions whose Fourier transform has a negative part and the potential diverges or has a cusp at 
zero. We also show that in the ground state configurations of the penetrable sphere model particles 
are superimposed on the sites of a close-packed lattice. 
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1 Introduction 



The present paper is a continuation of my earlier work on ground states of bounded Fourier-transformable 
pair potentials [1, 2]. Bounded or intcgrable interactions appear in quantum physics for example in 
Bogoliubov's theory of the Bose gas, which is based on models expressed in terms of the Fourier transform 
of the pair potential. Such interactions play also a central role in classical soft matter physics [3]. Those 
studied in [1, 2] had a nonnegative Fourier transform of compact support. In 2007 Likos and his coworkers 
published a mean-field study of the case when the Fourier transform of the pair potential has a negative 
part [4]. If the system discussed in [1, 2] showed already very peculiar properties at high densities, 
according to these authors a partly negative Fourier transform induced an even more curious behavior, 
the particles having the tendency to form clusters on the sites of a lattice whose lattice constant would be 
determined by the (negative) minimum of the Fourier transform, and only the population of the clusters 
would increase with the density. The question then arises whether it is possible to prove rigorously this 
kind of behavior. When, in the fall of 2008, I presented preliminary results on this problem and on its 
quantum mechanical counterpart [5] at a Montreal meeting in mathematical physics, Valentin Zagrebnov 
kindly informed me that, without knowing it, I worked on the theory of coherent crystals, a subject 
promoted by Russian physicists a long time ago. 

In effect, maybe the first attempt at a theory of crystallization was due to Vlasov. In his paper [6] 
Vlasov studied the solutions of the today called Vlasov equation 

9f df 1 df , 

— = v • — H F- — (1.1) 

dt dr m dv 

for the one-particle distribution /(r,v, t). There is no collision term, but the equation contains a self- 
consistent force field 

F(r,t) = -^ J u(r-r') //(r',v,i)dvdr' (1.2) 

induced by a translation invariant pair potential u. Considered on the torus or in the entire space, 
this equation has spatially homogeneous static solutions /°(v). Vlasov investigated his equation after 
linearizing it about / . Besides many others, he asked the question whether there may exist spatially 
periodic static solutions. If yes, they could be associated with crystals. In his analysis, however, he 
arrived at the false conclusion that for the existence of such a solution J u(r) dr < must hold. Today we 
know that there can be no equilibrium (crystalline or other) phase under the effect of such an interaction: 
the system collapses, the grand-canonical partition function diverges in finite volumes [7]. The possibility 
of a static periodic solution for a (super)stable u is, nonetheless, there, and to see how such a solution 
emerges, we can follow Vlasov's reasoning almost up to the end. Given the velocity profile /°, one 
substitutes f(r,v,t) = /°(v) + (fi(r,v,t) into Eq. (1.1) and keeps only the terms linear in <p: 

dt dr m dv dr 

where 

$(r,t) = J 0(r,v,t)dv. (1.4) 

With the ansatz 

<A(r,v,i)=e-*- ik - r 9k (v), (1.5) 

one obtains 

flk(v)= y ^ / ^(v')dv', (1.6) 
m[k • v — ui) J 

where u is the Fourier transform of u. Integration over v yields 

l = ffi/Ji^dv (1.7) 
m J k • v — lj 
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which implicitly determines the dispersion relation w(k). With the choice /°(v) = h(v 2 ) this further 
simplifies, and for a static solution (u> = 0), k must satisfy the equation 



1 _ 2fi(k) 



J ^'(v 2 )dv. (1.8) 



m 

In three dimensions, for a bounded and sufficiently fast decaying h this becomes 



l = - 2 3f^d, (1.9) 



Because h > 0, we sec that a solution for k is possible only if u has a negative part! Making the 
substitutions /°(v) = h(v 2 ) and w — also in g k , we obtain an approximate static solution of Eq. (1.1) 
in the form 

/(r,v) = Mv 2 ) + ^(v 2 )^c k cosk-r, (1.10) 

with the sum running over the vectors k which solve Eq. (1.9), and the real c k (incorporating the 
constant (2/m)u(k) J <? k ) chosen so that / > 0. To obtain (1.10) we have supposed that w(x) = u(— x), 
which implies w(k) = u(— k) real and (1.9) holding simultaneously for ±k. If u is integrable, then u(k) 
is continuous and decays at infinity; thus, for any "natural" u the sum in (1.10) is finite. Moreover, 
u(0) > for a superstable interaction, therefore u(k) < can be only if |k| > ko > 0. Let now /° be the 
Maxwell distribution, i.e., 

h(x)=p(^) 3/2 e-^, (1.11) 



. 2tt , 

where p is the density and /3 is the inverse temperature. Then h'(x) = —(/3m/2)h(x), and (1.8) becomes 

-/3pu(k) = 1. (1.12) 
This equation has a solution if (3p > {f3p) c , where 

Wp)c = . * u . (1.13) 

mm{u(k)} 

If the minimum of u is nondegenerate, at the critical point the resulting distribution is spatially periodic. 
As (3p increases, the solutions for k shift away from the minimizer of u, and the solution of the linearized 
equation is, in general, almost periodic. With the Maxwell distribution (1.10) becomes 



/(r,v) = /°(v) 



, /3m ^ ' 
1 — c k cos k • r 



(1.14) 



In order to keep / nonnegative, as (i increases, the coefficients c k must scale as (i^ 1 . 

The right condition for the formation of "coherent crystals" appeared only much later, in papers by 
Kirzhnits and Nepomnyashchii [8, 9]. The term referred to (hypothetical) crystals of mobile particles as 
in Vlasov's theory, capable of ballistic motion which must be coherent if it preserves long-range order. 
The context was somewhat different, these authors were interested in the crystallization of a quantum 
liquid. The starting point was a Hamiltonian the ground state of which was determined in the Hartree 
approximation. This analysis showed that the lowest-energy solution of the Hartree equation can be 
periodic only if the Fourier transform of the pair potential is partially negative, and the periodicity 
is then given by the wave vector at which the Fourier transform is minimal. These authors, just as 
Vlasov, were fully aware of the extraordinary properties of coherent crystals, even more pronounced in 
the quantum than in the classical case, and described them almost in the same terms as Likos et al. [4]. 

The mean- field or density- functional approach of Likos et al, the partial differential equation method 
of Vlasov, and the Hartree approximation of Kirzhnits and Nepomnyashchii are effective one-particle 
theories, and all arrive at the same conclusion. One can have little doubt in their truth. The problem is 



3 



more difficult than the case treated in [1, 2], and a new idea is necessary to obtain some progress in the 
rigorous theory. The idea presented and exploited in this paper is that of an infinite-density ground state 
(IDGS). Working in a finite domain (on a torus here), the TV-particle ground states are arrangements of N 
points that minimize the interaction energy. To each iV-point subset one can assign a discrete measure, 
and an IDGS is the weak limit of a sequence of discrete measures associated with iV-particle ground 
states, when N goes to infinity. The ground state energy and energy per particle diverge in this limit 
but, if the interaction is integrable, the energy per pair is convergent. The limit is proportional to the 
best superstability constant, the largest number that can multiply N 2 /V in a lower bound on the energy 
of N particles in a cube of volume V. The two notions, that of an IDGS and of the best superstability 
constant, can be related via an energy functional written for normalized measures on the torus. It is 
shown that any IDGS is a minimizer of the energy functional, and the value of the minimum is the best 
superstability constant. The task is then to find these minimizers, because they can provide information 
on ground state configurations at high but finite densities. Most results on IDGS's will be obtained 
by writing the energy functional in Fourier representation; an exception is the penetrable sphere model 
presented at the end of the paper (for another example see Ref. [10]). Our main result is a confirmation 
of the conclusion of the earlier works [4, 6, 8, 9]: if the Fourier transform of the pair potential has a 
negative part, the distribution of particles in high-density ground states is nonuniform. In some cases, 
e.g. when the potential diverges at the origin or has a cusp there, this nonuniform distribution cannot 
be approached through Bravais lattices. Except for the penetrable sphere model, from this analysis we 
cannot predict the ground state configurations at high densities with the precision obtained in Refs. [1, 2]. 
The lack of precise information is true also for interactions with a strictly positive Fourier transform, 
but we can at least assert that, in contrast with the former, the asymptotic distribution of particles 
as the density increases tends to be uniform. Since there is a continued interest in the Gaussian core 
model [11]- [16], even this weak result may be of some value. 

It is to be emphasized that the limit of infinite particle density is not at all unusual in classical physics. 
Whenever a continuum theory is applied to describe a system of classical pointlike particles, tacitly 
this limit is used. Theories of classical fluids, the rigorous van der Waals theory of liquid-vapor phase 
transition [17, 18, 19], the continuum theories of droplet formation [20, 21] and liquid-gas interfaces [22] 
are based on different free energy functionals defined on continuous mass densities. These theories can 
be derived through scaling limits in which the particle density tends to infinity while the mass of the 
particles and the interactions among them tend to zero in order to keep the mass and energy densities 
finite. In this work we do not have to scale the mass because it enters only the kinetic energy which 
vanishes in classical ground states (gravitating systems are out of the scope of this study), but we do 
scale the interaction: considering the energy per pair of particles corresponds to scaling the pair potential 
by dividing it with N as N tends to infinity in a fixed volume. The crucial difference compared with 
the continuum theories cited above is the extension of the energy functional to discrete distributions. 
Starting with them and by controlling their convergence in the limit of infinite density we can obtain 
information on ground state configurations at finite densities. 

The paper is organized as follows. Section 2 fixes the conditions on the interaction and the basic 
notations, and contains the definition of the best superstability constant together with some preliminary 
results on it. Throughout the paper we deal with both bounded and unbounded interactions. The results 
will be more complete for bounded interactions, and the reason of this appears already in Sections 2 
and 3. For bounded interactions the self-energy is finite, and the energy per pair with and without the 
self-energy converges to the best superstability constant as N goes to infinity from above and from below, 
respectively. This permits to prove that the convergence is uniform in the volume - a fact that we do 
not know for unbounded interactions. In Section 3 we still show that the infinite-density limit of the free 
energy per pair is independent of the temperature and is the same as that of the ground state energy. 
In Section 4 we define the infinite-density ground states as weak limits of Dirac combs associated with 
iV-particle ground states, introduce the energy functional and prove two theorems, one for bounded and 
another for unbounded interactions. They show that IDGS's minimize the energy functional, and the 
minimum is the best superstability constant. From the Fourier representation of the energy functional, 
introduced in Section 4.3, we derive some stability conditions in Section 5. Section 6 contains a complete 
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description of the stationary points of the energy functional. Their knowledge is important because there 
can be minimizers among them. The truly challenging problem is to find the ground state configurations 
(GSC) in infinite space. The definition, based on local stability, is recalled in Section 7. Here we prove 
two lemmas, the first relating periodic configurations and minimizers of the energy density, the second 
establishing the relation between periodic GSC's in infinite space and GCS's on tori. The analysis of 
the problem in infinite space is continued in Section 8. We introduce the notion of an IDGS in infinite 
space and formulate as a conjecture a commutative diagram relating four objects: GSC and IDGS in 
finite domains and in infinite space. The main result of this section is Proposition 8.2 which gives a 
new expression of the best superstability constant. Sections 9-12 present applications. In Section 9 we 
give a rather complete description of IDGS's in the case when the Fourier transform of the interaction is 
nonncgative. Section 10 contains the main theorem proving the asymptotic uniformity or non-uniformity 
of GSC's in infinite space for the case of a strictly positive or partially negative Fourier transform, 
respectively. In Section 11 we discuss two special classes of interactions with a partly negative Fourier 
transform. In the first one for any nonzero k the sum of the Fourier transform over integer multiples of 
k is nonnegative. In the second case the Fourier transform is ultimately positive and slowly decaying, so 
that the pair potential either diverges at zero or has a cusp there. In both cases we show that no Bravais 
lattice can be an IDGS, and high-density GSC's are different from Bravais lattices. On the contrary, 
bounded pair potentials that are flat or nesting at the origin or have a dominantly negative Fourier 
transform at large wave vectors may prefer the accumulation of particles on the sites of a lattice which is 
independent of the density. In Section 12 we demonstrate this property on the so-called penetrable sphere 
model in which the interaction is a repulsive square core potential. The paper ends with an Appendix. 



2 Best superstability constant 

Consider the problem of the ground state of a system of classical identical particles confined in a fixed 
bounded domain A of volume V and interacting via a translation-invariant pair interaction u(x — y) = 
u(y — x). Assume u to be integrable, bounded outside the origin, strongly tempered, superstable and 
lower semicontinuous (see below) . In the simplest case A is a d-dimensional cube of side length L taken 
with periodic boundary conditions, and u is also made periodic. The usual way to achieve this is to 
replace w(x) by 

u A (x) = u(x + Ln). (2.1) 

This series is absolutely convergent for any L > and any x (outside the set L1 A when u(0) = oo) if u 
is strongly tempered. The general definition of strong temperedness is 

Kx)|<qx|- d -" (2.2) 

for |x| > r , with some r ,C, r\ > 0. If u has no hard core, r = can be taken. Condition (2.2) also 
guearantees that u\ converges to u pointwise as L tends to infinity. A ground state of N particles in A 
is any TV-point configuration (r)jv = (ri,. . . ,rjv) minimizing the A^-particle interaction energy 

U A (r) N =^u A {r i -r j ). (2.3) 

i<j 

The minimum energy will be denoted by Eq(N). Because of the shift-invariance of u, any translate of a 
ground state is a ground state. To make sure that (2.3) can be minimized for all N, we will ask u to be 
lower semicontinuous [23]. 

A way to understand the behavior of soft matter at high density is to study the limit of infinite 
density. Since the particles have no hard core, this limit can be given a meaningful definition. The limit 
p = N/V — > oo will be realized by letting N diverge in a fixed volume. To characterize a ground state, a 
suitable scaling of Ua is necessary. If u is bounded then 

<r) N = -V—rU^T)* (2.4) 
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remains bounded in this limit whatever the sequence of Appoint configurations is. Moreover, 

e N = mm e(r) N = J E (N) (2.5) 

(r)jvCA Jy [Pi — I) 

is bounded as N goes to infinity even if u is unbounded but is integrablc and bounded from below. Indeed, 
divide A into N cubes of volume V/N and choose (r) N to be the centers of the cubes. Then the second 
sum in 

^w^^EE^^- 1 "^ (2 - 6) 

i j^i 

is a Riemann-sum, therefore 

e[v) N = X - y«(x)dx + o(l) (JV^oo) (2.7) 

is an upper bound to ejv- We conclude that the right quantity to look at is e[v)^ that, with a slight 
abuse, we call the pair energy. 

The search for the ground state in the limit of N going to infinity is closely related to finding the 
best super stability constant for u. The notion of superstability was introduced by Ruellc [24, 7]. Stability 
means the existence of a constant B such that for any N the ground state energy is bounded below by 
—BN; superstability means that the ground state energy density (energy per volume) increases with the 
particle density at least quadratically. If u is integrable and superstable, to leading order in the density 
the increase is quadratic. 

Definition 2.1 The best superstability constant for a superstable interaction u is the supremum of the 
positive numbers C such that for any large enough cube A of volume V 

U A (r) N > CN 2 /V (2.8) 
holds for every N above a possibly A-dependent value and every (r)jv C A. 
Lemma 2.1 The best superstability constant for u is 

C[u] =liminfC A [ti] (2.9) 

where the limit is taken over cubes of increasing volume V and 

C A [u] = liminf e N . (2.10) 

N— >oo 

Proof. Take any C < C[u]. From the definition of liminf it follows that Ca[u] > C for A large enough, 
say A D Ac- It also follows that there exists some N\ such that 

^ e N > C if A D A c and A^ > N A (2.11) 

which is (2.8). On the other hand, if C > C[u] then there exist arbitrarily large domains A such that 
C > Ca [u] , and an infinite sequence Nj of positive integers (which may depend on A) such that 

lim e Nj - C A [u\. (2.12) 

Thus, 

e N] < C (2.13) 
for j large enough, which is the negation of (2.8). □ 
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Remarks, (i) In the definition of superstability [7] there is an additional term —BN on the right-hand 
side of (2.8). We can add such a term with any B > without violating the inequality, it even allows for 
the extension of the inequality to every N, cf. Eq. (2.31) below. Obviously, no choice of B permits to 
increase the coefficient of iV 2 and attain a superstability constant larger than C[u] defined by (2.9). (ii) 
Later on, A will be a general parallelepiped. The periodized interaction can be defined, and the analogue 
of Lemma 2.1 can be proven in the following form. Let A be any nondegenerate parallelepiped, and for 
s > consider sA = {sx : x e A }. Then the best superstability constant is liminfs^oo C s a M- The 
limit actually exists, 

C[u] = lim CaM (2.14) 

A— s-oo 

if A tends to infinity in the Fisher sense [7] , and is the same as that one obtains with the use of u instead 
of the periodized u\. Based on the strong tempcredness of the interaction, a proof in analogy with the 
proof of existence of the thermodynamic limit of the free energy [7] could be done, but in this paper 
(2.14) is considered as a hypothesis. 

An integrable interaction u has a bounded continuous Fourier transform decaying at infinity [23] . We 
will denote it by v. Thus, 

v(k) = J u(r)e- ik r dr. (2.15) 
If X)keA* l w (k)l < 00 then X^keA* w (k)e lk r is a continuous function and 

ua « = v E ^ k ) eik ' r ( 2 - 16 ) 

keA* 

almost everywhere; if u is continuous, equality holds everywhere. Here A* = (2-7r/L)Z d if A is a cube of 
side length L. Superstability implies that v(0) > and also tt(0) > and ua(0) > 0. 

Proposition 2.1 C A [u] < v(0)/2 and thus C[u] < v(0)/2. 

Proof. This follows from Eq. (2.7). An alternative proof is obtained by computing the average of the 
potential energy, 

^l UM) „-_mm^i. ( , 17) 

The minimum is smaller than the average, 

^ eA , = £o(w) ,»^i) (2 , 8) 

and hence 

C A [u]<v{0)/2. □ (2.19) 

The following simple observation helps to see the possibility of superposition of particles in some 
ground state configurations. Given (t)m, for any integer n > 1 we define the nM-point configuration 
(r)^ by n times repeating (t)m, 

T mM+ j=Tj, m = 0,...,n-l, j = l,...,M. (2.20) 

Then, for any bounded complex- valued function / on K d , 

1 nM M 
ij—l 1 

When u is bounded, we apply this identity to / = Vma/2. Introducing 

V , \ N-l . . Vu A (0) 

= 2N^ ^ Uh{Vl Vj) = e(r)N + ~m (2 ' 22) 
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and 

N ~ 1 Vu A (0) 

= ( r ^CA W(r)jV = " + ' (2 ' 23) 

and recalling that m\n for integers m,n denotes that m is a divisor of n, we obtain the following. 
Lemma 2.2 

w(r)Jr = w(r) M (2.24) 

and, hence, 

u Nl > lo N2 > wjv 3 > • • • if JVi|JV2|JV 3 |---. (2.25) 
Proof. Equation (2.24) repeats (2.21). Let (r)^ be a A^-particle ground state. Then 

to Nj+1 < u{r)%+ l/Ns = uj(v) Nj = w Nj . □ (2.26) 
Proposition 2.2 If u is bounded, 

CaM = mfwiv. (2.27) 

Proof. 

. N -e N = ±(^fl-e N ) (2.28) 
goes to zero as N increases, so the statement is 

liminf cjjv = inf cjjv- (2.29) 

N^oo N 

Now infjv con < liminfA^oo cu^; supposing a strict inequality, because of (2.25) there would be a sequence 
Nj tending to infinity such that 

lon- < WjVi < lim infw tv (2.30) 

J A->oo 

contradicting the definition of liminf. □ 

A direct consequence of Equation (2.27) is that 

C/A(r)AT>-^JV + C A My (2.31) 
for all N. Another consequence is as follows. 

Corollary 2.1 If uj(y)m — Ca[m], then (r)^ is a nM-particle ground state for all n>l. 

The simplest realization of this situation is provided by the penetrable sphere model and is presented in 
Section 12. For a less trivial example see Ref. [10]. 

3 Infinite-density limit of the ground-state energy and the free 
energy per pair 

Proposition 3.1 The sequence is convergent. 
Proof. We give two different proofs. 

(i) The first proof works for u bounded. In this case wjj — 6 a tends to zero with increasing TV, therefore 
the convergence of ca is equivalent to 

lim wjv = inf cja = C\\u\. (3.1) 
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Now Ca[w] is the smallest accumulation point of the sequence u>n, and suppose there is another one, ujq. 
Let S = |(wo — C\[u]) and define two subsequences, K rn and L n via the inequalities u>K m < Ca[w] +6 and 
co Ln > cjq — S; thus, to Ln — cj Km > 5 for all m, n. Both subsequences are infinite and, because of Eq. (2.25), 
IK\ e {K m } for all integers I > 1. Therefore, for any n there is an m such that < L n — K m < K\ and, 
thus, 

E (L n ) < E (K m ) + l -K Y {K x - 1 + £r ro )|MU • (3.2) 

The right member of this inequality is an upper bound to the energy of a L„-particle configuration that 
one obtains from a if m -particle ground state by adding L n — K m particles. Dividing by L 2 n jV we find 



Wi„ - U} K „ 



OiL- 1 ) (3.3) 



contradicting ujl„ — ^K m > S. 

(ii) The second proof works even if u is unbounded. It is based on the monotonic increase of the ground 
state energy per pair, 

EojN+1) Ep(N) 

N(N + 1) ~ N(N-iy [ ' 

see also Kiessling [25]. If u is integrable then ejy — VE n (N)/N(N — 1) is bounded above by v(0)/2, cf. 
Eq. (2.18), so its limit as N goes to infinity exists, 

lim ejv = supe^v = Ca[u]. (3.5) 
As to Eq. (3.4), let N > 2 and R be any sequence of TV + 1 points of A (repetition allowed). Then 

U A (R)= Yl «A(x-y) = ^ E E u A(x-y) = ^ T ^ U A (S), (3.6) 

(x,y)Cfl SCR,\S\=N (x,y)cS SCR,\S\=N 

because for any two-point subsequence (x, y) C R we can choose N — 1 different TV-point subsequences 
S C R containing (x, y). Since the sum over S has N + 1 terms, 

min LT A (J2) > m in y A (S) (3.7) 

i?CA,|ii|=Ar+l TV - 1 ScA,|S|=iV 

which is just (3.4). □ 

We have nowhere used the periodicity of the interaction, therefore Proposition 3.1 applies also to 
u instead of u\. It holds true also for unstable interactions, when Ca[u] < 0. If u is bounded then 
Propositions 2.2 and 3.1 yield 

supe^v = C\[u] = inf uj n . (3.8) 

Corollary 3.1 If u is bounded then the convergence of £[pVJ and uj^ pV j to C\[u] as p tends to infinity 
is uniform in V. 

Proof. 

< max{w LpyJ - Ca[u],C\[u] - £y P v\} < uy pV \ - ey P v\ < < (3.9) 

ip p 

if p and V are large enough. Here we used Eq. (2.28) and superstability for the third and strong 
temperedness for the fourth inequality. □ 

The uniform convergence of £[ p yj to Ca[u] probably holds true also if u diverges at the origin. If A is a 
cube of side L and R is any TV-point configuration then there is an x G Asuchthatdist(x,i?) > L/(2iV 1 /' i ). 
Using this fact, for a radial u one can easily prove that 



-i ( 1 
ejv < ejv+i < ew + cp u I - 



(3.10) 



9 



if A is large enough. Here p = N/V and c = 1 if w(x) is a monotonic function of |x| close to zero. The 
term p~ 1 u(p~ 1 / d /2) is the analogue of p~ 1 u(0). Because u is integrablc, it tends to zero as p goes to 
infinity. However, the right member of (3.10) may not be an upper bound to C A [u] and, therefore, it 
cannot play the role of ojn in the bounded case. 

It is natural to ask how the free energy per pair behaves as the density tends to infinity. Let 

Za.n = I ^ UA(r)N d ( r V = e-^W. (3.11) 
For the free energy per pair, /at(/3), we have 

f < f (8)- VFn ^ <"(°) ( *m 

e N < fsifi) = N{N _ 1} < — (3.12) 

where the upper bound results from Jensen's inequality. We are interested in the limit of /jv(/3) as j3 or 
N goes to infinity. The usual definition of the partition function is V N /N\ times the expression (3.11). 
However, computing —(V/f3N(N— 1)) \yiZa.n with the modified definition or with (3.11) yields the same 
result for both limits. At positive temperatures, as N goes to infinity, /jv(/3) — £/v remains positive only 
if the entropy is negative and of order N 2 . This means that the level set 

{(r) N e A N : U A (r) N < E (N) + N 2 /(3} 

should be of Lebesgue measure <~ e~ cAr2 with some c > 0. Under some natural conditions on the 
interaction we shall find the opposite result, an infinite-density limit of the free energy per pair that 
agrees with the limit of ejv, i.e., Ca[u]. 

Theorem 3.1 Let u satisfy one of the following conditions. 

(i) u is bounded and 

|«(x) - «(y)| < (|u(x)| + |«(y)|) (3-13) 
with some r > and < a < 1. 

(ii) There exist constants u > 0, r > 0, c > 0, < ( < d and < a < 1 such that 



«(x) > mo ( o ) if |xj < r (3.14) 

and 



Then 



K x) - u(y)\ < c(M*)\ + Hy)\) ( min J, xU ^ ro} ) • (3-15) 
lim f N (/3) = e N , Hm f N (/3) = lim e N = CM\ all j3 > 0. (3.16) 



Proof, (i) Using (3.13), 



where 



|u A (x) - u A (y)| < 2||u A || [^f 1 ) ( 3 - 17 ) 
||u A || =sup \u(x + Ln)\. (3.18) 



Let (r)jv be a ground state configuration [hence, U\(r) n = E (N)] and (r+<5r)jv = (ri+5ri, . . . ,rjv+5rjv) 
an arbitrary perturbation of it. Then 

U A (r + 5r) N -E Q (N) < \\u A \\ ^2max^l^ N 2 < r 1 (3.19) 
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if for all i 



Using Eqs. (3.11), (3.12), (3.19) and (3.20), 



Cd 



1 > e -f>N(N-l)[f N (P)-e N ]/V > 

~~ ~ e 



(r e p , N /2) d ^ 



V 



(3.21) 



where Cd is a dimension-dependent constant. Taking the logarithm, dividing by /3N(N ~ 1)/V and letting 
either /3 or TV go to infinity we find Eq. (3.16). 

(ii) Let (r)jv be a ground state configuration and let = — r^. Because of Eq. (2.18) and the lower 
boundedness of u, + Ln) < uiN 2 holds with some u\ > for any n g Z d . By decreasing, if 

necessary, r , (3.14) is valid also with u\ replacing uq and implies 

|r y + Ln\ > r N~^. (3.22) 
Let x n = Yij + Srij + Ln and y n = r,y + Ln where Jr^ = Sri — Srj, and choose 

isr ^ r i(zj^y <3 - 23) 

for every i, where U2 = 3citi. If U2PN 4 > 1, this inequality implies 

2 



|*r y I <£i\T A. (3.24) 



Then 

7*n 2 

min{|x n |,|y n |,r }> yiV"3=< (3.25) 

and by (3.15), 

|u(x n ) ~ u(y,)| < c(Kx„)l + | U (y n )|) ( 2 ^ Tii \ ) . (3.26) 

\r a N ~J 

When passing to ua we have to estimate ^2 n |u(x n )| and ^ n |u(y n )|- It is easily seen that only a single 
term can be large, the sum of the rest is of order 1 as N goes to infinity. The possible large term is 
bounded by u\N 2 . Summation over i,j brings in another factor A^ 2 , and for N or (3 large enough we end 
up with a generous upper bound 

U a (t + 8t) n -Eq(N) < u 2 ^maxl^i^) 7V 4+ A < p~K (3.27) 

V 1 r a J 

The rest of the proof is as in the bounded case. □ 

Remarks, (i) For the validity of (3.16) the continuity of u was essential. We shall see that for 
the penetrable sphere model at least the first of Eqs. (3.16) fails. Condition (3.13) is a sort of strong 
Holder-continuity (because w(x) decays as |x| grows) which guearantees the ordinary Holder-continuity 
of u\. In the case of a free boundary condition ordinary Holder-continuity would suffice, (ii) There is 
a gap between interactions that are bounded or diverge algebraically at the origin. Interactions with 
a logarithmic divergence, occurring in some cases between colloidal particles [26], allow particles to be 
much closer to each other and, thus, the ground state configurations to be strongly inhomogeneous. One 
cannot exclude that instead of (3.20) one should impose \Srij\ = 0(e~ cN ). Then, one could still prove 
the first of Eqs. (3.16) but not the second one, because of an entropy oc —N 2 at positive temperatures. 
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4 Infinite-density ground state and the energy functional in fi- 
nite volume 



4.1 Bounded interactions 

In this section we show that it is possible to obtain Ca[u] and a good approximation of high-density 
ground states via the minimization of an energy functional. For bounded interactions we can write w(r) jy 
in the form 

= W /. UA<x ~ y ) d ^( r )«W d ^( r )«( y ) ( 4 - 1 ) 

where 



A 2 

N 

V 



is a measure on A of total weight V (henceforth, a normalized measure); 5 X is the Dirac delta at x. We call 
M(r)jv the measure associated with (r)^. The discrete measures of the form (4.2) - sums of Dirac deltas 
with equal weights, sometimes called Dirac combs - are dense among the normalized Borel measures 
relative to the weak topology, cf. Lemma A.l. Therefore, when N goes to infinity and Rn is a A-particle 
configuration, the sequence /ir n of associated measures can converge weakly to any normalized measure. 
We shall use this property to introduce the notion of an infinite-density ground state. Let M\ denote 
the set of normalized Borel measures on A. A sequence [i n £ Ma converges vaguely to a p e A^a if for 
every real continuous function / of compact support 

lim Hn{f) = M/), (4-3) 

where fj,(f) := J f d[i. The convergence is in the weak sense if (4.3) holds for every bounded continuous 
/. Because A is compact, every continuous function is bounded and of compact support, so the two types 
of convergences are the same. Henceforth /i„ — /i will denote that /i is the vague limit of the sequence 
yu n . For the reader's convenience, in the Appendix we show that any infinite sequence fi n € A^a has a 
weakly convergent subsequence (i.e., Ma is compact in the weak topology), and any /j, £ M A is the weak 
limit of a sequence of discrete measures of the form (4.2). 

Definition 4.1 /j, £ A^a is an infinite- density ground state (IDGS) of ua if there is a sequence ^iR m of 
measures associated with configurations R m such that \R m \ goes to infinity with m and 

/J, Rm ^(i, lim e{R m ) = C A [u]. (4.4) 

m— >oo 

Note that R m may not be a sequence of ground state configurations. 
Proposition 4.1 The set of IDGS 's is nonempty. 

Proof. Given any sequence {Rn}j° of A-particle ground states, e(i?7v) = ew tends to Ca[u], cf. Eq. (3.5). 
Moreover, the sequence [ir n of associated measures has at least one weak limit point which is, therefore, 
an IDGS. □ 

Motivated by Eq. (4.1), we define an energy functional on A4a by the equation 

If 1 

!m = 7^7 J^ 2 u a(x - y) d/i(x) d/x(y) = — \i * n(u A ). (4.5) 



Here we use the notations of Ref. [27]: for a complex function /, /(x) := /(— x), Jl(f) := A*(/). The 
measure V^ 1 ^*^ is the autocorrelation of \i. It is defined on A — A (= A in the case of periodic boundary 
conditions when A is a torus) by 

p *£(/)=/ /(z)d( M *M)(z) := / /(x-y)dM(x)d/i(y), (4.6) 

JA-A JA 2 
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where / is any continuous periodic function of period cell A. Thus, 

M*M= / #x- y d/z(x) djLt(y) = / (5 x+y d/i(x)d/x(y). (4.7) 

JA 2 J A 2 

For general results about the autocorrelation and its Fourier transform see [27, 28]. Equation (4.1) can 
be rewritten as uj(r)jy — I\fM( r ) N ]- is meaningful also for unbounded interactions, but it can be finite 
only if u is a continuous measure. For example, for the Lebesgue measure A, I[X] = v(0)/2. In I[fj] one 
can recognize the interaction energy common to all mean-field theories, either quantum as the Hartree 
approximation [8], or classical as the mean-field theory of fluids [4, 12]. The difference in our use of it 
is that, while in mean-field theories \i is tacitly supposed to be absolutely continuous, here no a priory 
assumption is made on it. Let 

I° = inf 7[4 (4.8) 

For bounded interactions 

I°< inf I\ji (r)N ]=Mu N = C A [u]. (4.9) 

-<v,(rjjv JV 

Before proving the opposite inequality in Theorem 4.1, let us note that the correspondence between a 
finite configuration and the associated measure (4.2) is not one-to-one. The configuration (r)^ 

defined by Eq. (2.20) is different from (t)m whose n-times repetition gives rise to it. This is in contrast 
with 

M(r)£ f =M(r) M > ( 4 - 10 ) 

implied by Eq. (2.24). Thus, any discrete measure is associated with an infinity of finite configurations. 
This observation leads to the counterpart of Corollary 2.1. 

Proposition 4.2 Suppose that Hm m is a minimizer of the energy functional. Then for any positive 
integer n, (r)^ is a nM-particle ground state of u\. 

Proof. Because of the identities (2.24) and (4.10), 

uj{v) n M = 7[ M(r)M ] = 1° < J[M(x) n „] - ^(x)„ M (4.11) 

and, hence, 

U A (r) n M < C/a(x)„m (4.12) 
for any nM-particle configuration (x)„m- n 

In this way, if a minimizer of I[fi] happens to be a discrete measure concentrated e.g. on the vertices of 
a lattice inside A, then the particles pile up on the lattice sites already in finite-density ground states. 
For a A-periodic bounded Borel function / and a u e M\, let 

iM = ^*m- (4.i3) 

7[/i] = I UK [lA in this notation. 

Lemma 4.1 If u is bounded then 7[/i] is continuous in the weak topology, i.e., I[fi n } — > I[lA if Hn — t 1 - 
Proof. Let first u be continuous. If \i n — \i then [i n * jl^ — ^ [i * Jj,, therefore 

if / is continuous on A. Applying this to / = u A the result follows. Suppose now that u is bounded. For 
any e > we can choose a real continuous / on A and an integer m (depending also on /) such that 
|| U A — f\\oo < £ and \If[fi] — 7/[u„]| < e for n > m. Then (Ima^] — / UA [Mn]| < 3e for n > m and, hence, 

lim \I ua [/j] - I UA [Hn} \ < 3e. (4.15) 

n— >oo 

This being true for every e > 0, the above limit actually vanishes. □ 
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Theorem 4.1 If u is bounded, then 

1° = CaM, (4.16) 
the infimum 1° is attained, and I\p] = 1° if and only if ji is an IDGS. 

Remark. Because of Eq. (2.27), Equation (4.16) states that the infimum of I[fi] over all [i <G A4a is the 

same as the infimum over discrete measures associated with finite-density ground states. This does not 

mean that the IDGS cannot be a continuous measure, see the forthcoming sections. 

Proof. Choose ii m £ Ma such that I[/x m ] tends to 1° as m goes to infinity (/i m may be the same for each 

m). Because the set of discrete measures (4.2) is dense in Ai\ and / is continuous, /i m can be chosen to 

be of the form (4.2); hence, I[/J, m ] > Ca[u]. If \i is any weak limit point of \i m then by the continuity of 

h 

C A M > 1° = I™ I\pT] = I[fj] > C A [u}. (4.17) 

This proves (4.16) and that /j, is an IDGS. Starting with an IDGS fi and choosing fi m — nn m where R m 
is the defining sequence (4.4), Eq. (4.17) shows that \l is a minimizer of /. □ 

Remark. We have found that weak limits of iV-particle ground states are minimizers of /, but have not 
proved that every minimizer = IDGS can be obtained as such a limit. 



4.2 Unbounded interactions 

We restrict the discussion to pair potentials such that w(x) — >• oo as x — >• and u is bounded otherwise. 
Definition 4.1 and Proposition 4.1 are valid in this case. Because the notion of an IDGS is independent of 
the energy functional, the divergence of on point measures does not allow to directly conclude that 
IDGS's are continuous measures. Below we present a proof which is independent of /. 

Proposition 4.3 Let u be integrable and u(x) — > oo as x — > 0. Then any IDGS is purely continuous. 

Proof. Take any x € supp^i. Let R m be the defining sequence of fi, cf. Definition 4.1. Because 
m a(0) = oo, all points of R m are distinct. Choose any e > and let C be an open ball of diameter e 
centered at x . By approximating the characteristic function of C with continuous functions one can see 
that 

lim f, Rm (C) = n{C). 

It follows that for m large enough fiR m (C) > \^{C) and, because /x_R m ({x}) = F/|i? m | for x 6 R nil 

\R m nC\>^^. (4.18) 



As a consequence, 

U A (R m n C) > — — 1 inf u A (x) 



W \ 2V J |x|< e 

and 



(4.19) 



8V 2 W\R„ 



inf ua(x)-^} (4.20, 

X|<£ Z 



where u A (x) = — min{uA(x), 0}. Because u is bounded below and strongly tempered, supu A is finite. 
Dividing by |i? m |(|i? m | — 1)/V and letting m tend to infinity, 

or 

< 4V(2CA[ U ] + ^up uX ) 

lnf|x|<e«A(x) 
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The upper bound vanishes as e goes to zero, proving that /x({xo}) = 0. □ 

The extension of Theorem 4.1 to unbounded interactions is not immediate because IDGS's are weak 
limits of discrete measures which make I[fj] diverge. We therefore introduce auxiliary functionals which 
can give finite values on point measures. For a positive K let 

Ik[/A = I «a(x - y) d/x(x) d/x(y) (4.23) 

W J{(x,y):« A (x-y)<K} 

and 

I[fj]= lim /kM=sup/kM. (4.24) 

If u is bounded, = I\p]. If u is unbounded, I[fj] < I[fj]. In contrast to /, / can be finite on point 
measures if u(x) tends to infinity when x — > but u is bounded outside the origin, in which case we also 
have 

HfA = ^7 lim / u A (x - y) d M (x) d/i(y) (4.25) 

where c?A(x,y) = min neZ d |x — y + Ln\. The usefulness of I is based on the easily verifiable fact that 
I[^m n ] — e(r)jv, and on 7[/x] = I[fj] (maybe infinite) for continuous measures. 

Lemma 4.2 is continuous then = < oo. 
Proof. Let 

9n = 7^(P x /i)({(x,y) e A x A : « A (x - y) > n}), (4.26) 
a monotone decreasing sequence. By the definition of the Lebesgue integral, 

I[H]= lim \l K [ l j\ + Kg K \=I[n]+ lim Kg Kl (4.27) 



K- 



K— ¥00 



implying that the limit exists (< oo). If jj, has a point part then 

2V lun^ g K = ( M x M ) ( f| {u A > n} J = ^ M ({x}) 2 > 0, (4.28) 

\n=l / xGA 

therefore J[/z] = oo, but is finite if e.g. fj, is concentrated on a finite number of points. Suppose that 
/j, is continuous, then lim gK = 0. If lim Kgx = 0, we have the claimed equality. If lim Kgx > c > 0, 
introduce 

7+ [ M ] = I K M - Jo M = M + /„- M • (4-29) 
Because ma is bounded below, the second term is finite. Now 

K-l K K' 

I K [tA > - 9n+l) = ^2(9n ~ 9 k) > ^{9n ~ 9k) (4.30) 

n— 1 n— 1 n— 1 

for any K' < K, therefore 

K' 

hm7+M>^ 5 „ (4.31) 

n=l 

for any K' . It follows that 

oo 

lim 7+M > ^5n = oo (4.32) 

71=1 

because <?„ > c/2n for n large enough. Thus, in this case = I[fi] = oo. □ 
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Lemma 4.3 If \i is an IDGS then < Ca[u]. 

Proof. Let R m be the defining sequence of \i. By Lemma 4.1 and Theorem 4.1, 

T K [l4= lim Ik[hrJ< lim I\HR m \= lim eLR m ) = C A [u]. (4.33) 
Taking the supremum over if gives the result. □ 
Lemma 4.4 C A [u] < /°. 
Proof. Define 

uf (x) = min{u A (x), X}. (4.34) 

For any AT there exists some -ftjv such that Rm is a iV-particle ground state of U\ if and only if it is also 
a TV-particle ground state of Uff for any K > Kn, and 

E*(N) = U«(R N ) = U A (R N ) = E (N). (4.35) 

For example, Kn = Eq(N) + (A^ 2 /2)supu^ is a possible choice. Indeed, < %, therefore Eff(N) < 
Eq(N). It follows that if K > Kn and Rn is a iV-particle ground state of cither U\ or Uff , then 

uf (x - y) = u A (x -y) < K 

for every pair {x, y} C Rn, otherwise 

U a (R n )>UK(R n )>E (N). 

This means, however, that Rn is a common ground state of Ua and fT^, and (4.35) holds true. With 
the above Kn and Rn, in an obvious notation, cf. (4.13), 

I°k n = C a [u Kn ] > = e N , (4.36) 

"a 

the three relations holding due to Eqs. (4.16), (3.5) and (4.35), respectively. Letting iV tend to infinity, 
or taking the supremum over N, 

1° = I° UA > lim I\ N > C A [u], (4.37) 

IV— >-00 "A 

where the first inequality follows from Eq. (4.34). □ 

Combining Proposition 4.3 and Lemmas 4.2, 4.3 and 4.4, for any IDGS fi we find 

7[ M ] = C K [u] = 1° = (4.38) 

and, therefore, the following result. 

Theorem 4.2 If u is integrable, bounded outside the origin, and u(x) — >• oo as x — > 0, then 1° = C\[u], 
the infimum is attained, and any IDGS is continuous and minimizes I. 

Remark to Sections 4-1 o,nd 4-2. In the proof of Theorems 4.1 and 4.2 we have not used that ma was a 
periodized interaction. Therefore, the theorems are valid in the case when A is any compact Lebesgue- 
measurable set and the interaction is the original, non-periodic u. 
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4.3 Fourier representation of the energy functional 

Most of the results on IDGS's will be found by writing in Fourier representation. This is obtained 
by substituting the expansion (2.16) into Eq. (4.5) and integrating term by term: 



M = \ 



v(0)+ J2 «( k Mk)| 5 
o^keA* 



1 jjL*Ji{v). (4.39) 



2V 



Here 

M(k) = l^e- kx d M (x), (4.40) 

so that |ju(k)| < /x(0) = 1, and 

1/1^= J2 \m)\ 2 5*- (4.41) 

keA* 

We do not need the equality (4.39) to hold for all fj, e A4\. We are interested in IDGS's, and I[fj] = 
Ca[u] < v(0)/2 for any of them. It suffices, therefore, that the Fourier expansion (4.39) converges to 
I[fj], whenever < v(0)/2. In the case v > 0, I[fi] < v(0)/2 holds (with equality) if and only if [i is 
an IDGS; later on, we will discuss this case in detail. Let 

i>+(k) = max{0,w(k)}, iT(k) = - min{0, t>(k)}. (4.42) 

If v~ 7^ 0, then the IDGS's are among those n making the right member of Eq. (4.39) absolutely convergent 
and satisfying 

£ v-(k)\m\ 2 - £ v+(k)\m\ 2 >0. (4.43) 

O^keA* O/keA* 

For a stable interaction the left member of this inequality cannot exceed v(0), cf. Proposition 5.1. The 
verity of (4.39) in this case can be seen as follows. Given a t > 0, consider the Gaussian 

G*(x) = (47rf)- d / 2 e-l x l 2 / 4t (4.44) 
and its periodization, called the periodic heat kernel, 

G A (x) = £ G*(x + Ln) = I £ e -*l k lV k - x , (4.45) 

n6Z d keA* 

where the right member is obtained by the Poisson summation formula (or the Fourier expansion of the 
periodic function on the left). If /i € M\, then G l A M\ as well. [/ * /i(x) := J /(x — y) d/x(y).] 

Moreover, G\ * ^ is absolutely continuous, and its Radon-Nikodym derivative 

0(x) = £ e- 4 l k l 2 /i(k)e lk - x (4.46) 

keA* 

is a real entire function of each component of x. These properties guarantee that the energy functional 
evaluated on G l A * /j, satisfies (4.39), 

i[g\*a = \y, <^y~ 2tM2 \m\ 2 - (4.47) 

keA* 

Because the left member is finite and the right member is absolutely convergent at t — 0, the regularity 
of Abel summability implies that Eq. (4.47) survives the t — > limit and yields (4.39). 
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The form (4.39) of I[fj] reveals an important property of the minimizers, hidden in the defining 
equation (4.5). For stable interactions < 1° < v(0)/2 but, because A* is a set of density ~ L d , for most 
H G M\ one finds of the order of some power of L. For example, 

i[vs x ] = \Y. v ^ = Vu ^°)/ 2 - ( 4 - 48 ) 

keA* 

Thus, one may expect that when A increases, the minimizers are among /x's such that for some volume- 
independent constant c and for any K > 

supp/iP|{k e A* : |k| < K } < cK d . (4.49) 

In Section 6 we shall meet some special sets of this kind, namely, additive subgroups of A*. These 
will be shown to support the Fourier transform of the stationary points of /. When v > or u(x) > 
const x e _Q ' x ' , at least the uniform local boundedness as A ^ oo of both the autocorrelation and its 
Fourier transform (4.41) is easily seen. We return to this question in Section 8. 



5 Stability conditions for pair potentials 

From Theorems 4.1 and 4.2 we can deduce stability conditions on u. If u is nonnegative then it is 
obviously stable; if u > and is integrable then its Fourier transform v is a positive definite function (a 
function of positive type) [29], therefore v(0) — \v(k)\ > for all k € K d . Due to Theorems 4.1 and 4.2, 
similar conditions can be obtained without assuming u > 0. 

Proposition 5.1 Let u be an integrable pair potential with u(—x) = u(x). Then 1° = I°(A) > for any 
cube A is sufficient and necessary for stability, and /°(A) > for all A is sufficient and necessary for 
super stability. Some particular necessary conditions for the stability of u are as follows, 
ft) 

2v(0) +v(k) > 0, allkeR d . (5.1) 
(ii) If u is radial and v(k) < for |k| > fc then 

v(0)+n d min v(k) > (5.2) 

|k|>fc 

where nd is the number of nearest neighbors in a d-dimensional close-packed lattice (nd — 6, 12, 24, 40, . . . 
for d = 2,3, 4, 5,...;. 

(Hi) If u is bounded, radial and u(r) < for |r| > ro then 

u(0) + n d min u(r) > 0. (5.3) 

r|>»"o 

Remark. Ruclle's example of a catastrophic potential ([7], Section 3.2.3) violates the inequality (5.3). 
Proof. Because for u bounded 1° = ini N ^ N w(r)jv, the condition 7° > is part of Proposition 3.2.2 of 
Ruelle [7]. The role of the sign of J° for u unbounded is also clear from 7° = Ca[u]. Now v(— k) = w(k) 
real, and assertions (i)-(ii) are nontrivial if v takes on negative values. 

(i) Suppose that w(k) < for some k ^ and consider the measure 

d^i k (x) = (1 + cos k • x) dx. (5.4) 

The form (4.39) of I[/j] shows (cf. also Proposition 9.1) that I[^k] = w (0)/2 + w(k)/4 which must be 
nonnegative if u is stable. 

(ii) Let min| k |> feo v(k) = v(k m ). Choose A and a Bravais lattice B (lattice, for mathematicians) such 
that the reciprocal lattice B* (dual lattice, defined by k • r e 27rZ for k e B* and r e B) is close- 
packed, B* C A* and the nearest neighbor distance of B* is |k m |. Then, for (r)w = B fl A, by a simple 
computation 

I\Mb] < ^M0)-n d Kk m )|], (5.5) 
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and the left member is nonnegative if u is stable. Here I[pb] = I[i^bda] which is independent of A 
provided that A is a period cell for B, see also Eq. (9.11). 

(iii) Condition (5.3) is the dual of (5.2), obtained from the Poisson summation formula 

«(r) = P{B) £>(k) = MB)I\mb] ■ (5-6) 

rGS keB' 

Here p(B) denotes the density of B. Let min| r |> ro u(r) = u(r m ). Choose B to be close-packed with 
nearest- neighbor distance |r m |. Then, in case of stability, 

° " = 2pJW) J u(r) " 2^B) [U{0) - n M*m)\] (5.7) 

which was to be shown. □ 



6 Stationary points of the energy functional 

The natural candidates for IDGS's are the stationary points of I\p\. To motivate the definition of 
stationary points, we first consider absolutely continuous normalized measures, 

d M (x) = / 2 (x)dx, (6.1) 

where / is a real function. We have to minimize I[p] with respect to / under the condition J A / 2 (x) dx = 
V. Via functional derivation of I[p] we obtain 



/(x) 



0. (6.2) 



/ w A (x- y)/ 2 (y)dy-c 

.J A 

Here c is a Lagrange multiplier. Another way to write this equation is 

"A(x-y)/ 2 (y)dy = c if /(x) ? 0. (6.3) 



I A 

Any absolutely continuous p with a density f 2 solving Eq. (6.2) is a stationary point of the energy 
functional. Because with such an / 

m = wL ua(x ~ y)/2(y)/2(x) dy dx = % (6 - 4) 

among the solutions of (6.2) one has to choose the one giving the smallest constant c. A solution of (6.2) 
is / (x) = 1 which defines the Lebesgue measure and yields c = v(0) and, hence, an already familiar 
upper bound 7° < v(0)/2, cf. Eqs. (2.19), (4.9). We shall see that the Lebesgue measure is the only 
absolutely continuous measure among the stationary points of the functional I. By generalizing (6.3) we 
arrive at the following definition. 

Definition 6.1 A /j, e Ma is a stationary point oflf, p G A4\ st (f), if 

(/ * /Lt)(x) = / /(x - y) d/i(y) = constant for x G supp^i. (6.5) 

J A 

fi G Ma is a stationary point (of I), p G *Ma s \ if it is a stationary point of If for every continuous 
function f. 
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Remarks, (i) One could ask / * fi to be constant only /x-almost everywhere, but the above choice is 
sufficient and more convenient for our purposes, (ii) The name is justified by the following property. Let 
fi £ MA St (f) and let v £ M\, v < \i. Then 

/ (/*/x)d/i= / (/*/i)di/, (6.6) 
ja J a 

and therefore 

J,[(l - e)n + ev\ - I f M = e 2 (I f [u] - I f \p\) (6.7) 

for any < e < 1. So fi may not minimize If, not even among the measures absolutely continuous with 
respect to it, but if we perturb /x with a v <C /x in the order of e, // changes only in the order of e 2 . 

MA St C A^a s *(/), and we shall give two examples showing that the first can be a proper subset of 
the second. MA St is nonempty, the Lebesgue measure is a stationary point. Intuitively, it appears to be 
clear that /x is a stationary point if it is uniform on its support, and supp fj, — x = {y — x|y £ supp fj,} is 
the same for all x £ supp /x. The theorem below provides the proof. 

For t £ A and a function g on A define T t g by 

T t <j(x) = <7(x-t), (6.8) 

and for zx £ Ma, It A 4 •= -Ma by 

T t n(g) = n(T- t g). (6.9) 

Then for Borel sets D c A, 

T tAt (£>)=7t M (l D ) = it(£>-t). (6.10) 
Lemma 6.1 Any translate of a stationary point of If is a stationary point of If . 

Proof. Suppose /x £ M.A St (f)- Let A — supp/x, and (/ * xx)(x) = c for x £ A. Then (/ * T t /x)(x) = c for 
xei + t and, because A + 1 = supply we conclude that T t x< € A / <A St (/). □ 

The set A = LT d is an additive group, A — A = A. A closed subset G of A is a closed additive 
subgroup if G — G C G. A Haar measure fi on G is a nonzero measure of support G which is invariant 
under shifts with elements of G; that is, T t xt = fi for every t £ G, \i is uniformly distributed on G. The 
normalized Haar measure is unique. 

Theorem 6.1 A fi £ Ma is a stationary point if and only if it is a translate of the normalized Haar 
measure on a closed additive subgroup of A. 

Examples. Some simple closed additive subgroups with their normalized Haar measures are 

{0} with VS , (6.11) 
A with \ d (6.12) 

and 

(^B^nA^xA^-^ (6.13) 

with 

^'"A X ^ E JnA(i)| E ^x^- (6-14) 

where n = 1, . . . , d, A„ and 6^ are the n dimensional Lebesgue and Dirac measures, respectively, A^™) is 
the projection of A to the first n coordinates and is a n-dimensional Bravais lattice commensurate 
with A("). Any Bravais lattice can appear as a stationary point if we replace A = LT d by a suitable 
parallelepiped. 

In the following lemma we give a full description of closed additive subgroups of A and the Haar 
measures on them. Related results can be found in Ref. [27]. Recall that A* = (2n/L)Z d and V = L d . 
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Lemma 6.2 Let G be a closed additive subgroup of A with normalized Haar measure fj, (i.e., [i(G) = V), 
and let 

G* = {k e A*|k • x e 2ttZ, aH x e G}. (6.15) 
T/ierc G* is an additive subgroup of A* , and Jl is a Haar measure (the counting measure) on G* , 

£(k) = ^ -5 k;K . (6.16) 

kgg* 

Conversely, if G* is an additive subgroup of A* with the counting measure fl on it, then 

G = {x e A|k • x e 2ttZ, all k e G*} (6.17) 

is a closed additive subgroup of A, and the measure fj, defined by ju as its Fourier transform is a Haar 
measure on G with /z(G) = V. 

Proof. The (closed) subgroup property is obvious in both cases. Let now G be a closed additive subgroup 
of A, cither given initially or constructed as in (6.17), G* the dual of G, and \i the normalized Haar 
measure on G. We show that ju is just (6.16). For any k e A* and y G G, 

m = \ J g e-*- x dM(x) = ^ J e" ik x d M (x + y) 

= 1 / e- 4k ( x - y ) d^(x) = e ik - y £(k). (6.18) 

Then either j2(k) = or k • y e 27rZ for all y e G, that is, k <G G*. In the second case j£(k) = 1, proving 
the claim. □ 

Thus, we have obtained a simple characterization of closed additive subgroups of the torus A as the 
dual sets, in the sense of Eq. (6.17), of intersections of linear subspaces of R d with the lattice A*. Indeed, 
each additive subgroup of A* is such an intersection. We now see that A is the only closed additive 
subgroup of positive Lebesgue measure, therefore the Lebesgue measure is the only absolutely continuous 
stationary point of /. 

Proof of Theorem 6.1. 1. Let [i be the normalized Haar measure on an additive subgroup G, extended 
to A with fi(A \ G) = 0. For arbitrary Borel sets Z?i C G and D 2 C A \ G, and any x e G, 

fi{D 1 -x)=fi(D 1 ), ^(D 2 - X )= f i(D 2 ) = 0. 

Hence, for any Borel set D C A and any x e G, ji(D — x) = ji(D) or, equivalcntly, 

(Id */u)(x) = n(D) for any x e G. (6.19) 

Because any continuous function is the limit of a sequence of finite linear combinations of characteristic 
functions, (6.5) holds for every continuous /; hence, (i € AiA St ■ 

2. Suppose that \i € A4\ st and let A = supp/z. Although the characteristic functions are discontinu- 
ous, Eq. (6.5) extends to them through limits. For example, the characteristic functions of open sets are 
limits of monotone increasing sequences of continuous functions. For an arbitrary Borel set D C A, 

Li(D + x) = \ D for any xe A (6.20) 

is then obtained by using the upper regularity of jj,. 

(i) Suppose that e A. Then /i(Z) + x) = /i(£>) for all D C A Borel sets and all x e A; in particular, 
for D = A - x, 

H(A - x) = fj,(A) = n(A n A - x) = V (6.21) 

for all x S A. Since A — x is closed together with A, necessarily A — x = A for all x <G A and, hence, 
A — A = A. Thus, A = supp \i is an additive subgroup of A and /i is a shift-invariant measure on it. 
(ii) If ^ A, choose any t e —A and replace A by A + t and fi by T t /z. Because T t fi € A4a s \ the 
argument (i) yields that T t M is the normalized Haar measure on the additive group A + t, and \i is a 
translate of this measure. □ 
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Corollary 6.1 Let /j,, v e M\ be two measures with coinciding supports. If both are stationary points 
then [i = v. 

Proof, supp^i = suppf is an additive group and \i and v are normalized Haar measures on it. But the 

normalized Haar measure is unique, whence /i = v. □ 

Remarks, (i) Here are two examples showing that M.A St ^ A^a^ («a)- 

1. A two-dimensional example is the Dirac comb fin concentrated on the honeycomb lattice. It is not a 
stationary point because, if Xi and X2 are nearest neighbors, supp fin — ~*-2 = — supp/Ujf +xi. However, 
u\(— x) = «a(x) implies fin € M\ st (u\). 

2. If u is radial and A = [— L/2, L/2] d , u\ has a cubic symmetry. Then for any < a < L/4 

is stationary for ma, but not for functions not having a cubic symmetry. 

(ii) If the periodic boundary condition is replaced by a free one, the set of stationary points changes 
completely The reason is that now u itself, and not the periodized ma, appears in the energy functional; 
in fact, A may not be a parallelepiped. The additive groups are replaced by others. For a radial u the 
stationary points of I u are uniform distributions on 1, 2,..., d-1 dimensional spheres, on the vertices of 
regular polyhcdra, in 3 dimensions also on the 60 vertices of the soccer ball (fullerene), and others. The 
supports arc now sets that are invariant under transformations forming finite point groups or groups 
of rotation, and the stationary points correspond to Haar measures on those groups. As A increases, 
the dependence of the ground states on the boundary condition should continuously disappear in the 
bulk. For example, on parallelepipeds there should be a continuous transition between IDGS's for free 
and periodic boundary conditions. Such a transition cannot take place on stationary points for both 
types of boundary conditions. Intuition suggests that stationary points have better chance with periodic 
boundary conditions. 

(iii) If u is integrable but u(x) — > oo as x — > then I\p] = oo for any [i having a point part. However, 
depending on the rate of divergence of m(x), I[fj] can be finite in some continuous stationary points, and 
these are to be considered as possible IDGS's. If m(x) ~ |x|~ d+< » near the origin with < ( < d then 
m is integrable in d — n dimensions for n < £, so the energy functional is finite on the corresponding 
measures (6.14). Some interactions not allowing stationary points of I to be IDGS's will be presented in 
Theorem 11.2. If < ( < 1, can be finite only if the Hausdorff dimension of supp/z is greater than 
d— £ > d— 1. According to our characterization given in Lemma 6.2, the only stationary point providing 
a finite energy for such an interaction is the Lebesgue measure. In Section 9 we will show, however, that 
the Lebesgue measure does not minimize I if the Fourier transform of u is partly negative. Thus, for 
m e L 1 (W i ) diverging at zero as |x| _d+ ^ with < ( < 1 and having a partly negative Fourier transform, 
no stationary point is an IDGS. Theorem 11.2 will cover this case. 



7 Ground state configurations in infinite space 

From Refs. [1, 2] let us recall the definition of a ground state configuration in R d . To begin with, a 
configuration is an equivalence class of point sequences, two sequences being equivalent if they differ 
only in a permutation of the entries. We must consider sequences instead of sets because particles can 
be placed in the same point if the interaction is bounded. For an infinite configuration X and a finite 
configuration R let 

U(R\X) = U (R) + I(R, X) (7.1) 

with 

U(R)= "(r-A I(R,X)= «( r - x )- (7-2) 

(r.r')C-R rG-R.xeX 
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Definition 7.1 Let u be a strongly tempered interaction. Given a real number ji, an infinite configuration 
X is a (grand canonical) ground state configuration for chemical potential /i (a fiGSC) if for any finite 
configuration Xf C X, U(Xf) is finite, the sum I(Xf,X\Xf) is absolutely convergent, and for any finite 
configuration R C R d , 

U(R\X \ X f ) - > U(X f \X \ X f ) - v\X f \. (7.3) 
X is a (canonical) ground state configuration (GSC) if (7.3) holds true for every R such that \R\ = \Xf\. 

Thus, a ground state configuration is a locally stable configuration, in the sense that no local modification 
can decrease its energy. Equation (7.3) must hold also when R or Xf is the empty set. Hence, if X is a 
AtGSCthen 

U(Xf\X\X f )-n\X f \<0 

for any finite Xf C X. If u is not strongly tempered, absolute convergence of I(Xf,X \ Xf) must be 
replaced by some weaker condition, e.g. Abel summability, see [2]. The above definition and the two 
lemmas below apply also to interactions that are not integrable at the origin or have a hard core. 

The existence of GSC's and ^iGSC's is by no means obvious. For certain interactions they can be 
constructed explicitly [1, 2, 10, 30, 31], for some others only their existence is proven [32, 33]. In this 
paper we assume that for the class of interactions considered, high-density GSC's do exist; at least, all 
statements about GSC's in infinite space are subject to this proviso. 

Definition 7.2 For x € R d , let C x denote the unit ball centered at x. An infinite configuration X is 
locally uniformly finite (l.u.f.), if there exists an integer m such that for any x e R d , 

\XnC x \<m. (7.4) 

Equivalently, for any compact K E R d there is an integer such that 

\XCi (K + x)| < m K , anyxeR d . (7.5) 

In [2] we proved the absence of metastability in the case of strongly tempered interactions, by showing 
that a necessary condition for a l.u.f. configuration to be a GSC is to minimize the energy density among 
l.u.f. configurations of the same density. Here we prove that this condition is also sufficient among 
periodic arrangements (while it is obviously insufficient among l.u.f. configurations). 
The general form of a periodic configuration is 

J 

x= l+JOB + y,), (7.6) 

where B = Zai + ■ • -+Zad, &i E R d are linearly independent, yj are not necessarily different c?-dimensional 
vectors, and in this paper we use 1+J to indicate that in the union coinciding points occur with repetition. 
The choice of B can be made unique by asking J to be minimal. B thus obtained is the maximal group 
of period vectors of X, and is denoted by B(X). 

The particle density and the energy per particle of an infinite configuration X are defined, respectively, 

as 

, v . .. \xnn\ u(xnn) 

p(X) = hm - — — — -, eJX) = hm — (7.7) 

provided that the limits exist. Here is a bounded Lebesgue-measurable domain and V(Ct) is its volume. 
The energy density (energy per volume) of X is e(X) = p(X)e p (X). If the particle density is kept fixed, 
we can work with any of e p or e. If X is a periodic configuration, then p(X) = \X n A|/V r (A) where A is 
any period-parallelepiped. If, moreover, u is strongly tempered, cf. Eq. (2.2), and f2 tends to R d in the 
Van Hove sense [7], then 

U(X ) 1 " .UjXo) I(Xq,X\Xq) 

epix) - T^r + pbi ^ ( °' m) - + — pbi — • ( } 
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Here X m (m = 0, 1, . . .) are non-overlapping translates of any finite Xo C X such that X = U^ =0 X m . 
The simplest choice is Xo = (yi, . . . ,yj), cf. (7.6). However, X n can be the content of an arbitrarily 
large period cell, and then comparison of (7.8) and (7.7) shows that 

I(X ,X\X )=o(\X \). (7.9) 

Since u is superstable, for some b > and p(X) > b/C[u] we have 

^<-^Mm (7 - 10) 

and, consequently, 

I(X ,X\X ) 

lim 1 " ' \ 0) = 0. 7.11 

The first part of the following lemma was proven by Sinai for a lattice gas [34] and by Radin in one 
dimensional continuous space [35]. 

Lemma 7.1 Let the interaction be strongly tempered, and let X be a periodic configuration, (i) If X 
minimizes the energy density among periodic configurations of the same density and of period vectors 
belonging to B(X), then X is a GSC. (ii) If X is a GSC then it minimizes the energy density among 
l.u.f. configurations of the same density (provided that the respective densities exist). 

Proof, (i) Suppose that X is not a GSC; then there exist Xf C X and R C R d finite sequences such that 
RHX = 0, \R\ = \X f \, and 

A - U(R\X \ X f ) - U{X f \X \ X f ) < 0. (7.12) 

We construct a periodic Y whose period vectors form a subgroup of those of X, such that p(Y) = p(X) 
and e p (Y) < e p (X). For an odd integer n, let 

X Q = Xr\A n , A n = i^Y^XiRi : - Xi < ^2~~} ' iy7 ' U " > 

This is the union of the content of n d unit cells of X. Define Y = X \ Xf U R. Choose n so large that 
Xf U R C A„. Let Y be the periodic extension of Y a , 

Y= (J (y„+t)= l+J (nB + y), (7.14) 

tenB yeYb 

to be compared with 

X= \J (X +t) = |+J (nB + x). (7.15) 

tenB xex 

It is straightforward to verify that 

A = U(Y \X\X )-U(X \X\X ). (7.16) 
Clearly, p(X) = p{Y) and from (7.8) we find 

\X \[e p {Y) - e p (X)] = A+ i [I(Y ,Y \ Y ) - I(Y , X \ X )] + l - [I(X ,X \ X ) - I(Y ,X \ X )} . (7.17) 

Now X \X f = Y \R implies 

/(Xo, X\X )- I(Y ,X \ X Q ) = I(X f ,X \ X ) - I(R, X \ X ), (7.18) 
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and, using I (A, B + t) = I(B, A - t), 



I(Y ,Y \ Y ) - I(Y , X\X ) = I(R, Y\Y )- I(X f ,Y \ Y ). (7.19) 

The four infinite sums in the right member of the last two equations are absolutely convergent and tend 
to zero as n tends to infinity, because the distances of R and Xf to X \ X n and to Y \ Y diverge with 
n. Thus, if n is large enough, e p (Y) < e p (X). 

(ii) Earlier it was shown (Proposition, Ref. [2]) that if X and Y are two l.u.f. configurations of equal 
density and e p (X) > e p (Y), then X is not a GSC. Since periodic configurations are l.u.f., the second part 
of the claim follows. □ 

For a general parallelepiped A spanned by the vectors = 1, . . . , d), let 

u A (r) = u ( r + t) (7-20) 

t£B A 

where 

B A = jx> i l i :ni,...,n d ezj. (7.21) 

Lemma 7.2 Let u be strongly tempered and let X be a periodic configuration, (i) If X is a GSC of u, 
then X fl A is a ground state of ii A on any period parallelepiped A of X . (ii) If X fl A is a ground state 
of u\ on any large enough period parallelepiped A of X , then X is a GSC of u. 

Proof, (i) If X is a GSC, then it minimizes e p among periodic configurations of density p(X). In 
particular, if R C A, \R\ = \X fl A| and Y is the periodic extension of R from A, then e p (X) < e p (Y). 
Comparison with Eq. (7.8) written for a periodic configuration Z of period cell A and X = Z n A shows 
that 

< iZnA):=E TzT^ = e ^ z) - 1 2 £ ( 7 - 22 ) 

Applying (7.22) with Z = X and Y, one finds 

e$(R) - e£(X n A) = e p (Y) - e p (X) > 0. (7.23) 

(ii) We prove that if Y is a periodic configuration, B(Y) C B(X) and p{Y) = p(X), then e p (Y) > e p (X). 
Let A be a period cell of Y (and, hence, of X) chosen so large that X n A is a ground state of ma. Then 
inequality (7.23) holds true with R = Y n A. □ 

8 Infinite-density ground state in infinite space 

Let M denote the family of (positive) Borel measures on R d . 

Definition 8.1 Let X n be a sequence of infinite configurations with existing density p n and energy per 
particle e p (X n ) which tend to infinity with n. Suppose that 

lim e p (X n )/p n = C[u] (8.1) 

n— ¥oo 

and there is some nonzero p, S M. such that 

VX n = P^ 5x ^ ( 8 ' 2 ) 

xex„ 

Then p is called an infinite- density ground state in M d . 
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In contrast with IDGS's in finite volumes, the existence of IDGS's in infinite space is not a priori 
guearanteed. The ultimate goal would be to prove the following commutative diagram: 



Conjecture 8.1 



GSCmk GSCm 



p—>oo 



p^°° (8.3) 



IDGS in A ► IDGS in R d 

A^co 

If u is superstable and strongly tempered then high-density GSC's and IDGS's in infinite space exist. Any 
IDGS in R d can be obtained as the vague limit of both a sequence of periodic extensions of IDGS's in 
increasing volumes and a sequence of Dirac combs (8.2) associated with GSC's of increasing density. For 
any high-density GSC X, the associated measure px is the vague limit of a sequence \xy n , where Y n is 
the periodic extension of a ground state of u\ n in some parallelepiped A n , and p(Y n ) = p{X). 

If the conjecture holds true then information on GSC's in R d can be obtained by studying IDGS's in 
finite volume. In the case when v > and of compact support the results of Refs. [1, 2] and those of 
Section 9 below prove the conjecture but nothing new can be obtained about GSC's. If v > but its 
support is noncompact, the existence of IDGS's in infinite space still follows from Ca[u] = C[u] = w(0)/2; 
the Lebesgue measure is one of the IDGS's (the only one if v > 0) and Section 9 gives a fairly complete 
description of them, cf. Corollary 9.1. Concerning the general case, below we present some clarification 
and partial result. Let A be any parallelepiped centered at zero, 



A = | ^ xA t : -1/2 < Xi < 1/2, * = 1, . . . , d j 



(8.4) 



where L are linearly independent vectors. If n is a positive integer, then nA = {nx : x e A} is the union 
of n d translates of A: 

nA=|J(A + t,) (8.5) 

where 



1 d 

tj = ^ y^mjjla, m,ji,mj2, ■ ■ ■ ,m,jd G {— n + 1, —n + 3, . . . , n - 3, n - 1}. (8.6) 



i=l 



Definition 8.2 Let n be a positive integer. The periodic extension E n from Ma to M n A is a map 
assigning 

n d 

£„M = ^lA+t,7t^ (8.7) 

to [i e Ma; that, is, 

E nl i{p)=n{D-t j ), j = l,...,n d (8.8) 
for Borel sets D C A + tj. The periodic extension E from Ma to M is a map assigning 

E f i = U+tlt/i (8.9) 
teB A 

to fi e Ma- 

Lemma 8.1 Ep is the vague limit of E2 n +it l as n tends to infinity. 
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Proof. Let / e C c (R d ). Comparing (8.9) with (8.7) one obtains 

Efi(f) = lim E 2n+ll i{f). □ (8.10) 

n— >oo 

Recall that A* = {J2 n 4 l- : m e Z, i = 1, . . . , d}, where 1- • lj = 2n5ij. Clearly, 

A* C (nA)* = n~ 1 A*. (8.11) 
Lemma 8.2 If n is a positive integer and k e A* then 

£rtk) = { ,gW (8.12) 
v ' y s(k)ju(k) 2/n «s even v ' 

w/iere s(k) e {— 1, 1}. J/ke (nA)* \ A* then E n p,(k) = 0. Furthermore, 

Ejj,(k) = l A .(k)At(k), k e M d . (8.13) 

Proof. Using the definition of i?„/i and (4.40), 

^ = L e ^ EM = E/ A ^ k - (x+t ^(dx) = ^x>-^ ( k). (8.i4) 

The average of_the complex units is 1, ±1, or 0, depending on the case. Taking the limit n — > oo gives 
the result for £/i(k). □ 

Below we use the notations 7 A [/z] and 7° for the energy functional (4.5) and its infimum (4.8), and 
£A(r)jv, £a,at, w A (r)jv and wa.jv for (2.4), (2.5), (2.22) and (2.23), respectively, to indicate the dependence 
on A. Expanding the sum that defines u A , 

Ia[(4 = \ I «(r) d7; 7A (r) = \e 1a {u) (8.15) 

where 7a = p * p/V for p e .M A . 

Proposition 8.1 For any positive integer n, C n m A [u](m = 1,2,...) is a monotone decreasing sequence. 

Proof. Let A m = n m A. From (8.15) or from Lemma 8.2 and the form (4.39) of the energy functional, it 
follows that for any p € M.\ m 

7 Am+1 [E n p] = 7 Am [p]. (8.16) 

By Theorems 4.1 and 4.2 then 

C Am+1 M=7 Am+i <7° m =C Am H. □ (8.17) 

For a stable interaction C A [u] > for all A, therefore the limit of C„™ A [w] as m tends to infinity exists 
and is bounded below by the best superstability constant C[u]. As noted earlier, the limit must, in fact, 
be C[u\. Now we prove a different form of C[u\. 

Proposition 8.2 If u is bounded then 

C[u] = lim - inf {eJX) : X periodic, p{X) = p} . (8.18) 
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Remark. Here we have not supposed that there are periodic GSCs for arbitrarily high densities. If this is 
the case, the infimum in (8.18) is attained, and C[u] = lim„_ ) . 00 e p (X n )/p(X n ) where X n is any sequence 
of periodic GSCs with p(X n ) tending to infinity. 

Proof, (i) First we find a sequence X n of periodic configurations such that p(X n ) diverges and e p (X n ) / p(X n ) 
tends to C[u]. Let A„ be any sequence of parallelepipeds such that Ca„[m] converges to C[u]. For each 
n let R nm (m = 1, 2, . . .) be a sequence of ground state configurations in A„, \R nm \ going to infinity with 
m. According to Equation (3.5), linim^oo e^ n (R nm ) — C\ n [u\. Therefore, one can choose a subsequence 
R nmn = R n such that limeA„(i?n) = limu>\ n (R n ) = C[u]. Let X n denote the periodic extension of R n . 
Thus, R n = X n n A„ and \R n \/V(A n ) = p{X n ) — > oo. From Eq. (7.22), for any periodic configuration X 
of period parallelepiped A 



eJX) |XnA|-l 1 , . <x u(0) 

1m = V?V £A{X n A) + W) u(r) = WA(X n A) - ^pW) • (8J9) 



Applying this formula to X n and A„, we obtain lim^oo e p (X n ) / p(X n ) = C[u\. 

(ii) Suppose there is a sequence Y n of periodic configurations with respective period cells A„ such that 
p(Y n ) tends to infinity and 

lim ?ljM = c Q <C[u}. 

We can always choose the sequence A„ so that it tends to infinity in Fisher's sense. We may also suppose 
that Y n n A„ is a ground state of u\ n , otherwise we can replace Y n by the periodic extension of an 
\Y n n A„|-particle ground state of u\ n , cf. Eq. (7.23). For any e > there is a N such that for any 

n > N, 



£p(Y n ) _ q 



p(Y n ) 



£ 

^3' 



e An (Y n nX)- ep(Yn) 



£ £ 



P(Y n ) 

The second and third inequalities use Eq. (8.19) and Corollary 3.1, respectively. By assumption, 

£ A„ (Y n H A) = e An ,|r n nA|, 

therefore 

\C An [u] - C \ < £ tfn>N. 
Thus, Ca„[u] converges to Co < C[u], but this contradicts the hypothesis (2.14). □ 

Remarks, (i) Without the hypothesis (2.14), from (2.9) we still find that equation (8.18) holds true if 
linip^oo is replaced by liminfp^oo. (ii) The first equality of Equation (8.19) applies also to unbounded 
interactions. Equation (8.18) is presumably valid in this case as well, but in the absence of a proof of 
Corollary 3.1 for unbounded interactions we cannot prove it. 

For the sequence X n used in the proof of Proposition 8.2, Eq. (8.1) holds true. However, the proof of 
Eq. (8.2) is missing. We may suppose that A„ = (2l n + l)Qo, where Qo is the unit cube centered at and 
the integers l n tend to infinity with n. X n can be selected so that \X n n Qo\ — max t6R d \ (X n + 1) n Qo\, 
then px„(Qo) > V{Qo) = 1 because px n (^n) = V(A n ). The task would be to prove that px n {Qo) is a 
bounded sequence. If this holds, one could find a vaguely convergent subsequence of px„ tending to a 
nonzero measure which is by definition an IDGS in R d . The same IDGS in R d could also be obtained as 
the vague limit of a sequence of IDGS's p n e Ma„ chosen such that p n (Qo) > ^(Qo)- 

Let G be a locally compact Abelian group (in this paper G is a torus or G — R d ). A measure p, on G 
is called shift-bounded if for any compact A C G there is a constant ola such that for all a <E G, 

p(A + a)<a A . (8.20) 

Suppose again that A„ is a sequence of increasing parallelepipeds such that C\ n [u] converges to C[u]. 
For any n let p n be an IDGS in A„. Hence, with the notation j n = V(An)^ 1 p n * p^, 

i lim 7„(uaJ = \ lim E-y n (u) = ^ lim Ej n (v) = \ lim %(v) = C[u). (8.21) 
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E"f n is periodic and, thus, shift-bounded. £/y„ is a positive measure which is the Fourier transform 
of a measure, therefore it is also shift-bounded, see [27] (Proposition 4.9) or [28] (Proposition 3.3). A 
necessary condition for (8.2) to hold is that these sequences are uniformly shift-bounded. In a special 
case the proof is given in the following lemma. 

Lemma 8.3 If u{x) > ce~ Q l x l 2 for some c,a > 0, then both sequences E^ n and £/y„ are uniformly 
shift-bounded. 

Proof. Let A be a unit cube. If 5 (k) = Ke-^/ Aa with K = ^ e - q2 / 4a dqj * then 1 A *9> U- / 
denoting the inverse Fourier transform of /, 

l^(k)| 2 = El n (A + a) < E^ n (l A+a * g) = E ln {l A+a g) < E ln (\l A+a \g) 

keA*nA+a 



< 



for n sufficiently large (/? = e ( 27r )d|y/a) d / 2 )• The passage to any compact A is obvious, so £/y„ is uniformly 
shift-bounded. From (8.22) or with a similar independent argument one can show that for any Schwartz 
function /, E^ n (f )(n = 1, 2, . . .) is a bounded sequence. We use this fact to prove that Ej n is uniformly 
shift-bounded. Let again A be a unit cube and g(x) as before; then there is some constant c g such that 

^7„(^ + a)<S 7n (l A+a *.g)-^(T A+a g)<^;(|U|g)<^;(g)< Cg all n. □ (8.23) 

An example to this lemma is the Gaussian pair potential for which, as mentioned already, the unique 
IDGS is the Lebesgue measure. Other examples are pair potentials with a Gaussian lower bound and a 
partly negative Fourier transform - the latter condition holds, for example, if Au(0) = 0, cf. Eq. (11.18). 
Unfortunately, uniform shift-boundedness of the sequences Ej n and Ej n does not allow one to conclude 
that the sequence fj, n is uniformly shift-bounded. Unbounded but o(y/V(A n )) weight rearrangements of 
fi n do not affect uniform shift-boundedness of Ej n . 

9 Infinite-density ground states for v > 

The most complete information about IDGS's can be obtained when the Fourier transform v of the 
interaction is nonnegative. The minimum of I\p] is attained on any fi not contributing to the sum in 
Eq. (4.39), and its value is 1° = v(0)/2, thus, C A [u] = v(0)/2 = C[u] for every A. 

Proposition 9.1 The Lebesgue measure A is an IDGS in A if and only if v > 0, and then it yields 

I[X]=C[u]=v(0)/2. (9.1) 

lfv>0, the Lebesgue measure is the unique IDGS. 

Proof. For k e A*, 

A(k) - <5 k , , (9.2) 

and A is the only element of Ma with Fourier transform <5k,o- So A is an IDGS, and is the only one if v 
is strictly positive. It remains to prove that the Lebesgue measure is not an IDGS if v is partly negative. 
In this case consider 

dyu q (x) = (1 + cosq • x) dx (9.3) 
and choose q € A* such that v(q) = mhikeA* w(k) < 0. Then 

/^(k) - S K0 + \ (<5k, q + 5 k ,_ q ) (9.4) 
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and, recalling that v(— k) = v(k), 

/ .<, W = H|) + £M < !|) =/[A1 , (9 . 5) 

so A is not an IDGS. □ 

Let us make a small detour here to explain the choice of /i q . The energy functional is the diagonal 
part of a quadratic form on M\. The operator underlying it has nice properties. 

Proposition 9.2 Let u be an integrable even pair potential, u(— x) = u(x). On L 2 (A) define an integral 
operator A by setting 

(A0)(x) - f UA (x-y)0(y)dy. (9.6) 

J A 

Then A is a bounded self-adjoint operator with eigenvalues {v(k) : k <G A*} and eigenvectors V-'k(x) = 
exp(ik • x). 

Proof. If k e A* then 

u(k)= y*u(x)e- tk x dx = y u A (x)e- tk x dx. (9.7) 
Substituting = -0k into Eq. (9.6) gives A?/>k = v(k)V>k- c 

Thus, ^ q is an absolutely continuous measure composed of the lowest-lying real eigenvector of A that we 
must mix with the constant eigenvector in order to make the sum nonnegative. 

If v > and takes on the zero value, besides A many other IDGS's appear. We consider two specific 
cases. 

Proposition 9.3 Suppose that v > and v(k) = for |k| > K , where K < oo. Let A be a paral- 
lelepiped, B any Bravais lattice such that A is a period cell for B, and let qs* denote the nearest-neighbor 
distance of B* , the reciprocal lattice of B. If qs* > Ko then 



\B n A 

1 1 xGBnA 



is an IDGS. 



Remarks, (i) For the finite point configurations associated with /xsnA, the proposition partly repeats the 
results of [1]. It provides also an example to Proposition 4.2. (ii) yU^nA extends periodically to R d into 

fi B = p(B)- 1 J2 *x, (9.9) 

xGB 

and the energy functional shown below actually depends only on [ib- 
Proof. For fisnA given by (9.8) and for k e A*, 

MsnA(k) = l B .(k)= ^.k- (9.10) 

KGB* 



This implies 

J[ Ms ] (9.11) 



/[mbha] = ^ 



O^KGB* 



The shortest nonzero vector in B* has a length qs» , at and above which w(k) vanishes, therefore /[^bha] = 
v(0)/2 = I°. □ 
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The second special case is an interaction with a "soft mode", v(k) — for |k| = fco and v is positive 
otherwise. One may think of the example 

w(k)= Cl e- C2 l k l(|k|-fc ) 2 (9.12) 

with some positive constants ci,C2,fco, but the precise form of v plays no role. Choose A so that A* 
contains a vector q of length fc ; then there are at least two vectors, ±q of length k . The following 
proposition is an immediate consequence of the discussion above. The measure /z q appearing in it is an 
example for an IDGS that is not a stationary point of /. 

Proposition 9.4 Ifv(k) > and v(k) = for |k| = fc , the measure /i q of Eq. (9.3) with |q| = fc is an 
IDGS. 

A particularity of the case v > is that convex combinations of IDGS's are also IDGS's. 

Lemma 9.1 If v > then I\pt\ is a convex functional on Ma, that is, if fi, u £ Ma and < a, (3 < 1, 
a + f3 = 1 then 

I[afj, + [3v} < al[fj] +f3I[v}. (9.13) 

Proof. This immediately follows from the form (4.39) of I\pt\ and the convexity of \z\ 2 on the complex 
plane: 

\a Zl +(3z 2 \ 2 <(a\ Zl \+(3\z 2 \) 2 <a\ Zl \ 2 + l3\z 2 \ 2 . □ (9.14) 
Proposition 9.5 lfv>0, the set of IDGS's is convex. 
Proof. Let /x, v £ Ma be two IDGS's and < a, (3 < 1, a + (3 = 1. 

1° < I[an + <3v\ < al[fi] + 0I[v] = v(0)/2 = 1°, (9.15) 
therefore afi + (3v is an IDGS. □ 

Corollary 9.1 If v > then the periodic extension of any IDGS from any parallelepiped A is an IDGS 
in R d . IDGS's in R d form a convex set. If v > 0, the Lebesgue measure is the unique IDGS in M d . 

10 Uniformity vs nonuniformity of high-density ground states 

Proposition 9.1 has an immediate implication on the distribution of particles in high-density ground 
states in A. If fi is an IDGS and R m is a sequence of A^ m -particle ground states in A such that [iR m 
weakly converges to it, then 

\_^nD\ = KD) 

m^oo p(R m ) 

for any open D C A, where p(R m ) = N m /V. Equation (10.1) is a direct consequence of the definition of 
an IDGS. If v > then A is the unique IDGS. Thus, for any sequence Rn of A^-particle ground states 
Hr n converges weakly to the Lebesgue measure A, resulting 

hm l^^l = \{D). (10.2) 

On the other hand, if v is partly negative, the asymptotic distribution of particles in any ground state is 
necessarily inhomogeneous. We shall prove an analogous result for ground state configurations in infinite 
space. 

In Rcfs. [1, 2] we gave a presumably complete description of defect-free high-density GSC's of inter- 
actions with a nonnegative Fourier transform of compact support. If the Fourier transform is strictly 
positive, a prominent example of which is the Gaussian pair potential, a rigorous identification of GSC's 
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is still missing and, if the decorrelation conjecture of Torquato and Stillinger [13] is valid, the task is near 
to impossible in high dimensions. One is then reduced to some probabilistic description, and such an ap- 
proach may be helpful already for d < 8. According to another conjecture of Torquato and Stillinger [14], 
for the Gaussian potential at d < 8 the close-packed Bravais lattice and its reciprocal lattice should be 
the unique GSC at low and high densities, respectively. However, this conjecture was disproved by Cohn 
and Kumar [15], who found uniform close-packed periodic structures of lower energy at low densities 
in 5 and 7 dimensions. The result below shows the fundamental difference between interactions with a 
strictly positive or a partly negative Fourier transform. It may also be helpful in the numerical search 
for low-energy arrangements at high densities. 

Theorem 10.1 For arbitrary d > 1, let u be a bounded, integrable, strongly tempered pair potential. 
Suppose that there exists a sequence of periodic GSC's X n of respective densities p n tending to infinity, 
and let 

Write X n in the form (7.6), 

J„ 

X n =\\)(B n + yj ), (10.4) 

where B n = B(X n ) = X)»=i Za„;. 

(i) Let the Fourier transform v of the potential be strictly positive. Suppose that the sequence of lattice 
constants |a„j|(i = l,...,d;n = 1,2,...) is bounded. Then [ix n converges to the Lebesgue measure in 
distribution sense, that is, for any Schwartz function f : R d — > C, 

lim J /(x) d Mx „ (x) = J /(x) dx. (10.5) 

(ii) If v is partly negative, then (10.5) does not hold for any subsequence of [ix n - 

Remarks. The restriction to bounded interactions is for technical reasons. The result implies that if 
v is partly negative, any (vaguely or in distribution sense) convergent subsequence of Hx n tends to 
an IDGS in infinite space which is different from the Lebesgue measure. The IDGS's are expected to 
be periodic or almost periodic in this case. Accordingly, the distribution of particles in high-density 
GSC's shows a pattern corresponding to the infinite-density limit. On the other hand, if v > then the 
distribution of particles is asymptotically uniform. We cannot prove that the set of lattice constants is 
always bounded, but the opposite, the divergence of at least one of the d lattice constants, leaves only 
two possibilities. The first is that J n , the number of particles in the primitive cell (called complexity), 
tends to infinity. The second is that J„ is bounded, implying that X n falls into the union of lower than 
d-dimensional substructures of diverging separation and density. In the case of the Gaussian potential, in 
5 and 7 dimensions, a new numerical result [16] suggests a uniaxially anisotropic high-density GSC with 
a stronger compression along the distinguished axis than perpendicular to it. According to part (i) of 
the theorem, the other lattice constants cannot remain bounded if the anisotropy is to survive the limit 
of infinite density. 

Proof, (i) Let A„ be a period parallelepiped of B n and, hence, of X n . Then, by Lemma 7.2, X n n A„ is 
a ground state of u\ n . From Eqs. (2.28), (3.8) and (9.1), 



On the other hand, 

^|x„nA n |=/[^J = ^ + ^ J2 «(k)lMx„(k)| 2 , (10.7) 
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where 



Thus, 



E y ( k )l^„(k)| 2 < — ^ < — ^ 

0#keB , P» P" 



for n large enough; moreover, for any bounded set D 

E ipx„(k)i 2 < — 



2u(0) 



inf ke£) u(k) ' 



(10.8) 



(10.9) 



(10.10) 



Take any / e <S(M d ) (function of rapid decrease). Using p n = J n p(B n ), the Poisson summation formula 
yields 

J fdpx^p- 1 E /(x) = J /(x)dx+ E /( k )^n(k). (10.11) 



Because the set of period lengths of B* (n = 1,2,...) is separated from zero, the sum over k ^ tends 
to zero as n goes to infinity. Indeed, for m e Z d , let 



Q m = {k € R d : m t < k, < to, + 1 (i = 1, . . . , d)} \ {0}. 



(10.12) 



Then U m Q m = R d \ {0} and \B* n Q m | < / with some / independent of n and m. Since u is continuous 
and strictly positive, for any M > there exists an e > such that w(k) > e if k G U| m |<M<3m- Write 



E = E + E 

|m| <MQm) kG-B* n(U j m | ;> jvfQm) 



(10.13) 



Using (10.10), by Cauchy's inequality 



E 7(k)px„(k) 

keB*n(U| m |< M Q m ) 



< 



2w(0) 

PnS 



E 



kG-B* n(U| m |< M Q m ) 

which tends to zero as n — » oo. On the other hand, 



|/>)|2 < 2u(0)ZMi/1g o 

Pn£ 



E /(k)Mx„(k) 

kes*n(U| m | >M Q m ) 



<< E 

|m|>M 



kGQ, 



l/(k)|. 



(10.14) 



(10.15) 



The sum on the right-hand side is convergent and tends to zero as M tends to infinity, proving Eq. (10.5). 
(ii) Fix Q = [-L/2,L/2] d and define N n = [p n L d \. Let Q n = [-L n /2, L n /2] d where L d n = N n /p n , 
then < L d - L d n < p- 1 . Let q <= Q* = (2ir/L)Z d be the minimizer of v(k) in Q* , and let q„ = 
(L/L n )q € Q n = (2tt I L n )Z d . Because L n — > L, q„ tends to q as n increases. Let G AIq^ be 
a Appoint approximation of pq n , cf. Eq. (9.3), constructed according to Lemma A.l. Let Z n be the 
periodic extension of R n . Then p(Z n ) = p n , and by part (ii) of Lemma 7.1, 



iim e p( x ") < lim e P ( z ") _ "(°) Kg) 



Pn 



Pn 



(10.16) 



Applying Eq. (8.19) to A„, 
e p (A„) u(0) 



^ + ^= WA „(A„nA„) = — 1— / / «(x-y)d/ix B (y)d/ix B (x), (10.17) 

'n <^Pn 2»/(A n ) J An J R d 
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where A„ is any period parallelepiped of X n . Suppose that fix n tends to A on Schwartz functions. If 
u G <S(R d ), this would yield v(0)/2 instead of (10.16). If u is not a Schwartz function, replace it by u*G l , 
compute the limit of e p (X n )/p n and let t tend to zero. However, 

f (u*G*)(x-y)dy = t;(0) (10.18) 

JWL d 

for any t > 0; thus, again, we obtain a contradiction. □ 

11 Interactions without Bravais lattice ground states at high 
density 

For certain potentials having a partially negative Fourier transform we can complete the information 
about GSC's given in part (ii) of Theorem 10.1. Namely, we will show that for these potentials no 
periodic GSC X can be a singly or multiply occupied Bravais lattice if p{X) is large enough, that is, 
J > 1 must hold and all y 3 cannot coincide [cf. (7.6)]. This result is valid in any dimension. Our method 
is to exclude Bravais lattices first as IDGS's and then as GSC's on the torus and in infinite space at high 
but finite densities. We note that one-dimensional examples of interactions having a nondegenerate GSC 
with two points in the unit cell [30, 36] and others with no periodic GSC [37] existed already thirty years 
ago, see also Radin's review [38]. 

11.1 Compensation by higher harmonics 

In the following example the negative contribution to the energy coming from a wave vector k is com- 
pensated by the positive contribution coming from integer multiples of k. 

Theorem 11.1 Consider a bounded integrable pair potential having a partly negative Fourier transform 
v with the property that for any k ^ 

oo 

^v(nk)>0. (11.1) 

n=l 

Then no Bravais lattice can be an IDGS: if B is a Bravais lattice and the parallelepiped A is a period 
cell for B, then I[pBnh\ > 1° ■ Furthermore, at finite, high enough densities there is no Bravais lattice 
among the ground state configurations on tori and in infinite space, and no lattice tower of the form 

J 

B J = |+J B (11.2) 

i=i 

where B is a Bravais lattice and J > 1 can be a GSC if p = J p(B) is high enough. How large p must be 
depends on the details of the interaction. As an example, suppose that 

(i) there exist < k\ < fc 2 < 2k\ such that v(k) < for k\ < |k| < fc 2 and v(k) > otherwise, 

(ii) v(2k) > \v(k)\ for k\ < |k| < k 2 , and 

(Hi) v(k) = for |k| > 3q with some q such that k\ < q < k 2 . 
Then there is no Bravais lattice GSC if 

p>{3/2) d -\<l/K) d - 

Proof. Consider Equation (9.11). If K € B* then nK € B* as well. Summing over B* along lattice 
half-lines, the contribution of each partial sum is nonnegative, which shows that I[pbda] > w (0)/2- On 
the other hand, choosing any q such that v(q) < 0, 7° < I[pq] = v(0)/2 — |w(q)|/4, cf. Eq. (9.5), thus, 
p BnA is not an IDGS. The absence of Bravais lattice ground states at high but finite densities is true 
both in finite and infinite volume. To see this, suppose that A is a period cell of B. Then 

I[p B ., nA ] = co(B J n A) = u>(B nA)>^. (11.3) 
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On the other hand, if p is high enough then according to Lemma A.l and due to the continuity of 
one can choose a configuration R = (ri, . . . , r| B ./ nA |) C A such that 



|/M-/[ Mq ]|<Kq)|/8 



and therefore 



I[H R ] = u(R) < 



v(0) Kq)| 



< uj(B j n A) - Kq)|/8. 



(11.4) 



This shows that B J n A cannot be a ground state on A. In infinite space we compare the energies per 
particle e p (B J ) and e p (i? per ) where i? per is the periodic extension of R. Now B J and R per have the same 
density, and Eqs. (8.19) and (11.4) imply e p (i? pcr ) < e p (B J ). From Lemma 7.1 we conclude that B J is 
not a GSC. 

Consider the specific example. The idea is to find a periodic discrete point approximation X of the 
density 1 + cosq • x such that w(X n A) < v(0)/2. The computation is done by using the expression 



1 1 N 

= 5 E «(k) JjT, 



keA* 



n=l 



(11.5) 



see Eqs. (2.22) and (2.16). Select a q whose length q satisfies condition (iii). Choose L so that q — (2ir/L)l 
where I is an integer. Let ej, j = 1, . . . , d be the Cartesian unit vectors chosen in such a way that q = qei, 
and 

A= |X> e J ■•■0<x i <L). (11.6) 



For a positive integer l±, with a = L/l± define the lattice 

d 



B„ 



2tt 
— i 

q 



-iei + a^] njej : (m, . . . ,n d ) E Z d 



(11.7) 



J=2 



and, for any integer M > 2, the periodic configuration 



M 



X= \J(B q +y m ) 



where 



M q 

With this choice of L and a, the reciprocal lattice 



m— 1 



171 - 1 7r 

y m = ^7 ~ e i> m=l,...,M. 



(11.8) 



(11.9) 



2tt 



S* = <| mqH y^rijej : (m, . . . , n d ) e 



J=2 



(11.10) 



is a subset of A* = (2w/L)Zr. We show that for p > (3/2) (q/n) the energy per particle of X is lower 
than that of any B J of the same density. Now 



1 



IX n A| 



E 



|s 9 n A| 



M 

y e ik-rJ_y- ik 



rSB,nA 
M 



m— 1 



(11.11) 



rn— 1 
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if k e A*, therefore the sum (11.5) reduces to B*. The contribution of k = to to is v(0)/2. For k = ±q 
we have 



-T, 



M 2 



M ' 6 

m=l 



M 2(l-cos^) 1 (f) , M^OO 



and the function decreases monotonically with M to its limit at M — oo. For k = ±2q, exp(ik • y m ) runs 
over the Mth roots of unity and, therefore, X)m=i e 4k ' y "* = 0. Under the condition that all other vectors 
of B* have a length not smaller than 3q, 

mnA) = ^- v (n.13) 

The density of X is 

Mq 

P = - ; - . 11.14 

In one dimension £?* = gZ, so the condition of the validity of (11.13) is satisfied, and the smallest density 
to which (11.13) applies is q/w, obtained with M = 2. The density can be increased only by increasing 
M, therefore the pair energy (11.13) increases with the density, but for all p > q/ir 

w( *nA)<^-0) Kq)|, (11.15) 

where the upper bound is obtained by substituting (11.12) with its limit at M = oo. For d > 2 we can 
reach the smallest value for both lo(X n A) and p if we choose the smallest M, that is, M = 2, and the 
largest a compatible with the condition of validity of (11.13). For a the condition implies 2ir/a > 3g, so 
its largest allowed value is a — 2ir/3q. To summarize, with (8.19) we find that the energy per particle of 
any infinite-volume ground state configuration G of density p > (3/2) d_1 (g/7r) d satisfies the inequality 

e P (G) < Of- - ^ - ^ < ep(B .) _ EhM (d > 2) ( ii.i 6) 

and 

.^.M.^.iM <,_,, (1L17) 

if B is any Bravais lattice and J > 1 such that Jp{B) = p. □ 
11.2 Potentials with a cusp at zero 

If the integral / w(k)|k| 2 dk is absolutely convergent, then u is twice continuously differentiable at 0, and 

-Au(0) = {2it)- d j v(k)|k| 2 dk. (11.18) 

If 

lim / u(k)|k| 2 dk = oo, (11.19) 

K ^°° J\k\<K 

then u has a cusp at zero. We consider as a special case of a cusp if u(r) — > oo as r — > 0. The condition 
(11.20) below implies that u has a cusp at zero. 

Theorem 11.2 Let v~ be nonzero and of compact support. Suppose that there exist K, c > such that 

«(k)>-^ for\k\>K. (11.20) 



3G 



Let A be any parallelepiped. Then, no Bravais lattice commensurate with A can be an IDGS in A. In 

particular: 

(i) If d > 1 and 

<v(k)<-J— at\k\>K (11.21) 



| k |d+2 - V > - | k |™+2+ £ 

for some positive integer n < d, positive d and < e < 1, then IDGS's may be singular continuous 
measures of the form 

MO = ^d-n X MfiMnAW , (11.22) 

cf. Eq. (6. 14-), in which case high-density ground state configurations in A are strongly anisotropic, and 
may or may not be Bravais lattices, 
(ii) If d> 1 and 

v(k) > |^jg for \k\ > K, (11.23) 

then no IDGS is a stationary point of the energy functional. Furthermore, high-density ground state 
configurations are different from Bravais lattices and their towers (11.2) in A and, if u is bounded, also 
in infinite space. 

Remark. At the end of Section 6 we discussed the case when tt(x) ~ |x| ^ d +C w ith < ( < d near the 
origin. By a Tauberian argument w(k) ~ 1 1<: | at infinity. If also v~ ^ 0, then d > 4 and 3 < ( < d is 
covered by (11.21), and < C < min{d, 3} is covered by (11.23). If < C < 1 then Y,n=i v ( n ' k ) = 00 for 
all nonzero k, which is a special case also of Theorem 11.1. 
One- dimensional example to (ii). Let 

v(k) = ^e- a / fe2 + (1 - k 2 )e- k2 (11.24) 

For a > 3 this function takes on negative values on two symmetric bounded intervals and decays as 1/fc 2 
at infinity. Thus, u is bounded but has a cusp at zero, and 

E i£( fc )i 2 < 00 

must hold for an IDGS. Hence, any IDGS is absolutely continuous in A = [— L/2,L/2] and is different 
from the Lebesgue measure, because v~ ^ 0. High-density GSC's form a non-arithmetic sequence. 
Proof. Let us return to the tempered measure G A * \i and Equation (4.47). For any t > 0, 

^/[GW] = - E «(k)|k| 2 e - 2t l k l 2 |£(k)| 2 

= E ^(k)|k| 2 e- 2t l k l 2 |M(k)| 2 - E ,+ (k)|k| 2 e - 2t l k l 2 |£(k)| 2 , (11.25) 

keA* keA* 

because both sums are absolutely convergent. Now lim t _x) d/[G A * At]/di exists, it can be — oo for a 
general /i, but it must be finite nonnegative if n is an IDGS. Thus, for any IDGS [i in Ma, 

E ^ + (k)|k| 2 |M(k)| 2 < E ^(k)|k| 2 |/i(k)| 2 < E ^(k)|k| 2 < oo. (11.26) 

keA* keA* keA* 

Let B be a Bravais lattice commensurate with A and let fisnA be the associated measure, cf. Eq. (9.8). 
Then 

E « + (k)|k| 2 |^|(k)| 2 = E « + (k)|k| 2 = oo (11.27) 

keA* kes* 

because of (11.20), so [IbdA is not an IDGS. 
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Combining (11.26) with (11.20), we find that 



keA* 1 1 

must hold for any IDGS p. Thus, either j£ is nonvanishing on a lower than ci-dimensional subset of A* or 
it decays to zero. 

(i) In this case 

w+(k)|k| 2 <oo, (11.29) 

keA*:fc 1 = ... = fe d _„=0 

therefore a measure of the form (11.22) [which is a stationary point of /, c.f. Eq. (6.14)] cannot be 
excluded from being an IDGS. Indeed, 

]T «+(k)|k| 2 |M5(k)| 2 = ]T «+(O d _„,k')|kf <co. (11.30) 
keA* k'e(B<™')* 

Now /io can be obtained as the weak limit of a sequence of measures associated with Bravais lattices: 
the component p B in) is unchanged in this limit, while Xd- n is approximated by a sequence of discrete 
measures associated with Bravais lattices that fill in A( d ~") more and more densely. The approximating 
measures can, thus, be associated with more and more anisotropic Bravais lattices. 

(ii) Combining (11.26) with (11.23) one finds 



£ 

keA* 



(11.31) 



for any IDGS p. Apart from the Lebesgue measure, no stationary point of the energy functional satisfies 
Eq. (11.31). Moreover, in contrast with A and po, a p ^ A satisfying (11.31) cannot be obtained as the 
weak limit of measures associated with Bravais lattices. Indeed, for any Bravais lattice B, p^X(k) = 1 
on a sublattice of A*; therefore, if p ^ A was a weak limit of measures associated with Bravais lattices, 
then p(k) = 1 would be on lattice, and (11.31) would fail. It follows that there exists some pa such that 
if B is a Bravais lattice, A is a period cell for B and p(B J ) = J\B n A|/V > p\, then B J n A is not a 
ground state of u\. In one and two dimensions a pair potential can be bounded and satisfy condition 
(11.23). Then, due to the uniform convergence proven in Corollary 3.1, p\ can be chosen independent of 
A, and we can conclude that in infinite space B J is not a GSC of u. □ 



12 The penetrable sphere model 

So far we have only given one example, Proposition 9.3, illustrating Corollary 2.1 or Proposition 4.2. 
However, in that case GSC's of the form of lattice towers — particles superimposed on lattice sites — 
are submerged in a continuum of other GSC's. The truly interesting situation is when a lattice tower 
ground state occurs in a non-degenerate manner. Intuition to find such interactions may be based on 
both r- and k-space considerations. The k-space representation (4.39) of I[p] suggests that IDGS's are 
to be looked for among Dirac combs concentrated on Bravais lattices if v(k) < for |k| large enough. 
Choosing the density of the lattice B large enough, apart from pb(0) = 1, /As(k) will be nonzero only 
where i>(k) < 0, with its modulus equal to 1 on the lattice B* . Thus, one may expect that I[p] is 
minimized by some high-density ps, see also Eq. (5.5). If v has a unique minimum at k m , the simplest 
guess is that B* is close-packed with a nearest- neighbor distance (close to) k m . Thinking in r-space, 
the best candidates are potentials that are flat or even nesting at zero, that is, having Au(0) > 0. The 
interactions tt(x) ~ exp(— a|x|^) with j3 > 2, studied by Likos et al. [4], belong to this class. Intuitively, 
such an interaction may prefer the formation of lattice towers. The large negative part of v in such a 
case, seen on (11.18), supports this intuition. Below we discuss the penetrable sphere model which is the 
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simplest although somewhat pathological example. A one-dimensional family of pair potentials giving 
rise to superimposed particles in ground state configurations is presented in a separate publication [10]. 
The particles in the penetrable sphere model interact via the pair potential 

, . / u > for |x| < d Q ,.„., 
U W = ( for |x|>db. (12J) 

We define u to be lower semicontinuous. This potential is applied in soft matter physics to model a 
system of interpenetrating micelles in a solvent [39]- [42]. The interaction has a partly negative Fourier 
transform. In 3 dimensions it reads (k = |k|) 

Anuodl ( sin kd \ . . 

v(k) = — — cos kd . (12-2) 



(kd a ) 2 \ kd 

Note that v is not absolutely integrablc. The first and deepest minimum of t>(k) is at k = k m = 6.12/rfo- If 
this was to determine the periodicity of the ground state then the ground state would be the dual of an fee 
lattice of lattice constant k m , that is, a bec lattice of nearest-neighbor distance y/En/km = (7.70/6. 12)rfo, 
cf. Eq. (15) of Ref. [2]. Instead, we will see that the ground state is an fee lattice of nearest-neighbor 
distance do- 

We show that with a proper choice of A, the ground state configurations can be given for every 
particle number N. In d dimensions choose ai, . . . ,a<j such that |aj| = do and the vectors aj generate a 
<i-dimensional close-packed lattice, 

B = j^rija, : n u ...,n d e zj . (12.3) 

Choose A to be a period cell for B containing \B fl A| = M points of B and having a side length > 2do 
in every direction. The volume of A is 

V = M|det[ ai ...a d ]|, (12.4) 

so the volume of the unit cell of B is 

,--V/*-|*h....dl-{^£ (12.5) 

The property of B is that for the given density ps it has the largest nearest-neighbor distance or for 
the given nearest-neighbor distance d the largest density among Bravais lattices. We make also the 
hypothesis (trivial in two and proven in three dimensions) that B realizes the densest packing of hard 
spheres of diameter do. Note that with the choice of A we exclude the occurrence of other close-packed 
(e.g. hep) structures which can be ground states in suitable domains or in infinite space. 

If N < M, the ground state in A is continuously degenerate, any configuration with all inter-particle 
distances > do is a ground state. Below we describe the ground state configurations (GSC) for > M, 
modulo cyclic translations in A. 

Proposition 12.1 The IDGS is B n A, 

C A H=w(BnA). (12.6) 

For N finite, let n > 1 be an integer and let nM < N < (n + 1)M . Then there are ( N M nM ) GSC's in A. 
In each of them (n + 1)M — N points of B fl A are occupied by n particles and N — nM points of B fl A 
are occupied by n + 1 particles. Moreover, 



nV 

so that 



(n + l)M 



2N 



(12.7) 



n r l i- u ° u ° ~ / ^ d l/ 2 if d=2 , 1oa , 

GaM = hm ejv = = — x < . >- „ , „ (12.8) 

AL J tv^oo 2p B 2 \ dl/^2 if d=3. v ; 
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Proof. The proof goes by induction over TV. Because each side of A has a length > 2do, only a single 
term of ua can be nonvanishing, thus, u\(r) = u or 0. If N = M, B(lA has zero energy, therefore it is a 
GSC. Any perturbation of Bfl A creates at least one pair of particles of distance < do, therefore of energy 
uq. Thus, B is the unique GSC in A. Suppose we know the result up to N, where nM < N < (n + 
and want to prove it for N + 1 . Adding a particle to a ground state of N particles costs the less if it is 
placed on a point of B n A occupied by n particles. The increase in energy is nu ; at any other place it 
is at least twice as much. No relaxation of the configuration thus obtained can decrease the energy, so it 
is a ground state for N + 1 particles. The ground state energy for N particles is 

E (N) = (N- nM) n(n 2 +1) ^ + [(n + 1)M - iV]!^_ll Uo = nUo [ N _ ^M]. (12.9) 

Dividing by N(N - 1)/V yields e N . □ 

An unpleasant feature of the penetrable sphere model is the discontinuity of its energy at zero tem- 
perature when N > M. The values of Ua^n can only be integer multiples of u taken from the set 
{Eq(N), . . . , N(N - l)uo/2}. Therefore, 



1 
V~ N 



-I3[F N (P)-E„(N)] = J_ I e -l3[U A (r) N -E (N)] d ^ N 

A™ 



N(N-l)/2 

]T e -^ mu ° AdN ^ m) . (12.10) 



yN 
m=mo 

Here XdN denotes the (iA^-dimcnsional Lebesgue measure, 

A m = {(r) N € A N : U A (r) N = E (N) + mu } , (12.11) 

and 

too = min{m : \ d N(A m ) > 0}. (12.12) 

Now A is the set of translates of B n A increased by the finite number of possible assignments of 
particles to M sites, so that Xd(A ) > 0, but its higher than d-dimensional Lebesgue measures vanish. 
On the other hand, 

^n(A n{n ^)/2) = ^A d ( J B do/2 ) JV - 1 > 0, (12.13) 
where B do/2 = {x e R d : |x| < d Q /2}. Therefore < to < N(N - l)/2. From (12.10) it follows that 

c -0m o u o Arfjv(A m J e -i3[F N (l3)-E„(N)] < e -/3m «o _ (12.14) 

Taking the logarithm and dividing by — /3, 

m Q u < F N (f3) - E Q (N) < m u Q + i In [^/A^^J] . (12.15) 
Letting /? tend to infinity in (12.15) and in the mean energy 

W) ___ W)t ^__. (12 . 16) 
we find the following. 

Proposition 12.2 For the penetrable sphere model in a finite volume, 

lim F N (0) - lim E N (fi) = E (N) + m u , (12.17) 

0— >oc 0— >oo 

where too is given by Eq. (12.12). 

From the point of view of thermodynamics, it is, therefore, more adequate to consider E (N) + m u as 
the ground state energy. The limit of /jv(/3) as N goes to infinity depends on the limit of mo/N 2 and on 
the large N behavior of XdN{A mo ), and remains to be answered. 
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Note added 



After the submission of this paper I learned that a part of the tools I developed for studying ground 
state configurations at high densities existed already in abstract potential theory. The first result dates 
from 1923 and is due to Fekete [43]. In an article about the distribution of the roots of polynomials of 
integer coefficients, Fekete considered the quantity d n (S), the maximum of the geometric mean of the 
n(n — l)/2 distances between points in an n-point set, chosen from an infinite compact set S. He proved 
that d n (S) was monotone decreasing as n increased, and called the limit the transfinite diameter of S. If 
we take the pair potential w(x — y) = — In |x — y|, then the ground state energy per pair of n particles 
in S is — In d n (S). Thus, Fekete's result is the proof of the monotone increase of the ground state energy 
per pair [equation (3.4)] and of Proposition 3.1 in this special case. What I called the best superstability 
constant is minus the logarithm of the transfinite diameter. The first proof of the monotone increase 
of the ground state energy per pair for a general symmetric kernel u(x, y) was given by Choquet [44], 
who also proved a theorem corresponding to the present Theorem 4.2. The best superstability constant 
appears in [44] under the name of Fekete constant. The analogous theorem for a kernel which is bounded 
on the diagonal (Theorem 4.1 here) can be found in a recent paper by Farkas and Nagy [45]. Early results 
on abstract potential theory are summarized in Fuglede's paper [46]. 
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Appendix 

In this appendix we collect a few results on the convergence of measures on compact sets. Let A C M d 
be compact and let Ai\ be the set of normalized Borel measures on A. 

Lemma A.l Any ji S A^a is the weak limit of a sequence {hr n }^ = i of discrete measures of the form 

Proof. We may suppose that A is a cube, otherwise, we take a cube that covers A and extend fj, with zero 
value to a measure on the cube (cf. Lemma A. 4.) We may also suppose that A is a unit cube. Given 
/i G Ma (thus, /j,(A) = 1), divide A into a disjoint union of n d semi-open cubes Q™ of side length 1/n, 

where n = [log 2 (N + 1)J . In Q™ choose N 3 points {r"J^\ where N 3 = LN>(<9")J or Nj = \Nfi(Q])~\ in 
such a way that V . Nj = N. Thus, 



|MQ?) - N 3 /N\ < 1/N. 



(A.1) 



Define 



, n d Nj 



3 = 1 t=l 

where S r is the Dirac measure at r. For any continuous function /, let 



</>nJ = 



(A.2) 



(A.3) 



Because A is compact, / is uniformly continuous, so there is some 5 n which goes to zero as n increases, 
such that for each j 



1 N > 

^ -</>»„ 



<s n . 



(A.4) 
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It follows that 

n d 

\(f)n,j\- (A.5) 



The second term is 0(n d /N), so both terms tend to zero as iV goes to infinity. Thus, fiR N /i. □ 

Next, we show that given any infinite sequence of normalized measures on a compact set, one can 

select a subsequence converging weakly to a normalized measure. We start with the compact set being a 

unit cube A equipped with periodic boundary conditions (e?-dimensional torus). 

The Fourier transform (4.40) ofa^e Ma is a function of positive type [29] on A* = 2nZ d , meaning 

that for any integer l 7 any ki, . . . , k/ in A* and any complex numbers z\, . . . , z\ 

i 

M(k* - kj)ziZj > 0. (A.6) 

It follows among others that 

|/x(k)| < £(0) = 1. (A.7) 

Lemma A. 2 From any sequence fi n G Ma one can select a subsequence [i ni such that fi ni converges 
pointwise to the Fourier transform fl of a fi € Ma- We say that [i ni converges to \i in Fourier transform. 
Proof. Because A* is a countable set, we can use the diagonal process [23] to choose a subsequence fl ni 
which converges pointwise to some function ip on A*, 

lim ju ni (k) = ^(k), every k G A*. (A.8) 

i— >oc 

Fixing any integer I, ki, . . . , k; in A* and complex numbers z\,...,zi, writing down Eq. (A.6) for J2 ni 
and taking the limit i — > oo, we find 

i 

V-(ki - kj-)«i«j > 0. (A.9) 

Thus, ip is a function of positive type, ip(0) — 1, and by Bochner's theorem [29], tp = ft, the Fourier 
transform of a probability measure /x on A. □ 

Lemma A. 3 [i n converges to [i in Fourier transform if and only if it converges weakly to [i. 
Proof. If fi n — t 1 then 

lim/i„(/) = fj,(f) for /(x)=sink-x, cosk-x, (A. 10) 

therefore lim/i„(k) = /2(k) for all kg A*. 

Suppose now that fi, [i n G Ma, V- = lim£2„. Consider the trigonometric polynomials, 



P (x) - ake *' x ( A - n ) 



keA* 

where =/= for a finite number of k. Any continuous function / on A is the uniform limit of a sequence 
p m of the form (A. 11). Given any e > 0, we fix an m so that ||/ — p m ||oo < e, and choose N E such that 
for n> N £ 

HPm) - »n(Pm)\ < ]T K P( k ) ~ £™( k )l < ^ ( A - 12 ) 
keA* 

Because — p m )| < e and |/z n (p m — /)| < e, we obtain 

<3e if n>7V £ , (A.13) 

proving that \i n converges weakly to /U. □ 

The extension of the above result to the general case, when the compact domain may not be a cube 
and the boundary condition may not be periodic, is immediate. 
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Lemma A. 4 Let Aq C R d be a compact Borel set and let denote the set of normalized Borel 

measures on Kq. Given any infinite sequence \i n € M.k q , one can select a subsequence converging weakly 
to some fi G A4a - 

Proof. Take a cube A D Aq. Let p! n e A4a with supp^ C Ao be the extension of jj, n to A. Constructing 
a normalized limit measure [i 1 G A4\ according to Lemma A. 2, it is clear that supp// C A , therefore 
M = m'|a c € A4a - Because // is the weak limit of [i' n , cf. Lemma A. 3, [i is the weak limit of [i n . □ 
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